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CHAPTER  I.  INTRODUCTION 


1 .0  Background 

The  basic  operation  of  a  radar  is  the  transmission  and  reception 
of  electrical  energy.  The  received  signal  or  radar  return  is  composed 
of  target,  noise,  jammer  and/or  clutter  energy.  For  a  ground-based  air 
defense  radar,  the  target  is  an  aircraft,  missile,  etc.;  clutter  is 
ground,  trees,  rain  or  chaff;  and  jammers  are  electrical  energy  trans¬ 
mission  devices.  Whereas  proper  radar  design  will  reduce  the  effects 
of  clutter  and  jammers  while  enhancing  the  target,  a  signal  processor 
is  normally  required  to  provide  target  enhancement  while  rejecting 
interference.  Additional  interference,  i.e.,  thermal  noise  from  system 
electronic  components,  increases  the  total  interference  power  which  the 
processor  must  reduce. 

A  typical  quadrature  channel  radar  digital  signal  processor  is 
shown  in  Figure  1.  The  mixers  and  the  lowpass  filters  are  used  to 
translate  the  intermediate  frequency  bandpass  radar  signal  to  in-phase 
and  quadrature  channel  baseband  signals.  After  the  signals  are  digi¬ 
tized  by  the  analog-to-digital  converters,  the  clutter  power  is  reduced 
by  clutter  rejection  filters,  commonly  referred  to  as  moving  target 
indicators  (MTI).  Once  the  clutter  power  is  reduced  below  the  thermal 
ncise  level,  the  signal -to-noise  (and/or  jammer)  ratio  is  improved  by 
some  type  of  coherent  integrator,  for  instance,  a  fast  Fourier 
transformer.  With  the  clutter  rejected  and  the  signal -to-noise  ratio 
increased,  the  amplitude  of  each  range/doppl er  cell  is  extracted. 

These  detected  outputs  are  sent  to  a  decision  element. 
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The  function  of  the  decision  element  is  to  produce  an  output  or 
target  report  only  if  a  target  is  present,  i.e.,  a  detection.  If  no 
target  is  present  and  an  output  is  reported,  this  is  a  false  alarm 
which  the  decision  element  should  minimize.  Typically,  the  probability 
of  detection  should  be  greater  than  50  percent  while  the  false  alarm 
rate,  or  probability  of  false  alarm,  would  be  between  10'°  and  10  . 

Basically,  the  decision  element  compares  a  threshold  (which  is  a 
function  of  the  system  probability  of  false  alarm  requirement)  to  the 
detector  outputs.  A  target  report  is  issued  if  the  threshold  is 
exceeded. 

If  a  fixed  threshold  decision  element  is  used,  the  false  alarm 

rate  is  extremely  sensitive  to  small  changes  in  the  average  value  of 

the  energy  from  all  sources  of  interference.  This  sensitivity  is 

easily  seen  in  Figure  2,  If  the  threshold  is  set  for  a  probability  of 
-8 

false  alarm  of  10”  ,  an  increase  of  only  3  dB  in  total  interference 

4 

power  density  corresponds  to  a  10  increase  in  the  probability  of  false 
alarm.  This  increase  would  put  an  unreasonable  demand  on  the  radar 
data  processor.  Therefore,  an  adaptive  threshold  decision  element  is 
required  to  provide  acceptable  target  detectability  while  maintaining 
a  constant  false  alarm  rate  (CFAR). 

The  processing  principles  used  to  counteract  the  variations  in  the 
output  interference  level  are  referred  to  as  constant  false  alarm  rate 
(CFAR)  or  adaptive  detection  processing  techniques.  The  most  common 
approach  to  the  design  of  such  CFAR  processors  is  to  sample  the  back¬ 
ground  interference  in  the  time-and/or-frequency  domain  around  the 
current  range  and  doppler  cell,  then  utilize  the  samples  to  estimate 


the  unknown  statistical  parameters  of  the  interference.  This  estimate 
is  used  to  maintain  a  CFAR  by  control  of  the  threshold  level . 


INCREASE  IN  NOISE  POWER  DENSITY  IN  OB  FROM  DESIGN  VALUE 

Figure  2.  False  Alarm  Probability  for 
Fixed  Threshold  Detection 


1 .1  Purpose 

The  purpose  of  this  study  is  to  compare  the  performance  of  two 
commonly  known  CFAR  techniques,  the  cell  averaging  and  the  "greatest- 
of."  The  performance  comparison  will  be  based  on  detection  probabili¬ 
ties  and  false  alarm  probabilities  obtained  from  a  Monte  Carlo  simula¬ 
tion  of  the  two  techniques.  The  performance  of  the  processors  will  be 
determined  for  different  target  models,  clutter  environments,  detector 
laws,  wordlengths,  and  interfering  target  power  levels.  The  simulation 
will  be  utilized  for  verification  of  theoretical  performance  results 


5 

and  determination  of  performance  results  not  obtainable  by  present 
analytical  methods. 

1 .2  Content 

Chapter  II  reviews  the  theoretical  analysis  of  a  fixed  threshold 
processor.  Probability  density  functions  at  the  detector  output  are 
derived  for  noise  only  and  target  plus  noise  for  both  a  linear  and  a 
square  law  detector.  Two  target  models  are  used,  a  steady  or  nonfluc¬ 
tuating  target  and  a  Swerling  I  target.  The  probability  of  false  alarm 
for  each  detector  is  determined  by  using  the  noise  only  probability 
density  functions.  The  probability  of  detection  is  determined  for  each 
detector  and  a  steady  target,  and  for  the  square  law  detector  and  a 
Swerling  I  target. 

Chapter  III  contains  a  theoretical  analysis  of  three  commonly 
known  adaptive  threshold  or  CFAR  techniques.  The  CFAR  techniques  are 
cell  averaging,  "greatest-of and  log.  Each  technique  is  described 
and  various  performance  equations  are  derived  or  given.  Probability 
of  false  alarm  and  probability  of  detection  equations  are  derived  for 
the  cell  averaging  and  "greatest-of"  CFAR  methods. 

Chapter  IV  describes  the  Monte  Carlo  simulation  developed  for  a 
performance  comparison  of  the  cell  averaging  and  the  "greatest-of"  CFAR 
techniques.  Mathematical  models  of  the  targets,  noise  and  clutter  are 
given.  Implementation  of  the  two  CFAR  techniques  is  presented. 

Finally,  the  determination  of  the  probabilities  of  false  alarm  and 
probabilities  of  detection  for  both  CFAR  processors  is  discussed. 

Chapter  V  presents  a  performance  comparison  of  the  two  CFAR 
processors  based  on  the  false  alarm  and  detection  probabilities  obtained 
from  the  Monte  Carlo  simulation.  The  comparisons  include  different 
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detector  laws,  clutter  environments,  quantization  wordlengths,  and 
interfering  target  power  levels. 

Chapter  VI  contains  the  summary,  conclusions,  and  recommendations 
for  future  work. 


CHAPTER  II.  FIXED  THRESHOLD  PERFORMANCE  ANALYSIS 


2.0  Introduction 

The  radar  statistical  detection  problem  in  noise  is  one  of  choosing 
between  signal  and  noise  at  the  radar  processor  output  or  noise  alone; 
that  is,  when  the  processor  output  voltage  is  described  by  v(t),  one 
wants  to  test  per  range  cell  between  H0  (noise  alone)  or  (signal 
plus  noise)  as  follows: 

Hq:  v(t)  =  n(t) 

H1 :  v(t)  =  s(t)  +  n(t)  .  (2.1 ) 

The  fixed  threshold  analysis  assumes  that  the  decision  element  is 
preceded  by  a  prewhitening  or  clutter  rejection  filter,  such  as  an  MTI. 

The  decision  element  in  Figure  1  tests  the  processed  video  to 
determine  whether  a  signal  is  present  (H^)  or  not  present  (Hq).  For  a 
specified  voltage  level  or  threshold,  the  decision  element  reports  a 
target  if  the  amplitude  of  the  video  is  greater  than  the  threshold,  and 
reports  no  target  if  the  video  amplitude  is  less  than  the  threshold. 

It  is  possible  that  processed  video  which  contains  only  noise  can 
exceed  the  threshold  generating  a  false  target  report  or  false  alarm. 

By  increasing  the  threshold  the  number  of  false  alarms  diminishes. 
However,  the  chances  of  detecting  a  target  also  decrease.  Consequently 
the  threshold  setting  is  made  as  low  as  possible,  consistent  with  a 
tolerable  false  alarm  rate  with  which  the  system  can  operate. 

In  a  given  system  the  fixed  threshold  would  be  determined  by 
establishing  a  tolerable  false  alarm  rate  based  on  overall  system 


considerations.  Having  determined  the  threshold  setting,  the  proba¬ 
bility  of  detecting  a  desired  target  can  be  calculated. 
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( 
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This  chapter  will  discuss  the  probability  density  functions  (pdf) 
at  the  output  of  the  detector  for  both  a  square  law  and  a  linear 
detector  with  noise  only  and  signal  plus  noise  inputs.  The  signal  or 
target  models  used  were  a  nonfluctuating  or  steady  target  [1]  and  a 
Swerling  I  target  [2].  Expressions  are  given  for  the  probability 
of  false  alarm  and  the  probability  of  detection  associated  with  the 
pdf. 

2.1  Noise  Only 

This  section  gives  the  pdf  and  probability  of  false  alarm 
expressions  for  the  single  pulse  amplitude  detected  noise  only  cases. 
The  detectors  considered  are  the  linear  and  square  law  detectors. 

2.1.1  Linear  Detector 

A  linear  detector  extracts  the  envelope  of  the  video  and  is  given 
as 

z  =  Jx 2  +  x2  ,  (2.2) 

where  Xj  is  the  in-phase  video,  Xq  is  the  quadrature  video,  and  z  is 
the  detector  output. 

If  Xj  and  Xg  are  independent  zero  mean  Gaussian  random  variables 

2 

and  homogeneous,  i.e.,  they  have  the  same  variance,  a  ,  the  pdf  of 
z  is 

p(z )  =  exp[-z2/2a2]  .  (2.3) 

o 

This  pdf  is  the  well  known  Rayleigh  distribution. 
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Assuming  a  fixed  threshold,  Z.^,  the  probability  of  false  alarm  is 


PFA 


/co  /*  °°  no 

p(z)dz  =  J  K  exp[-z2/2o2]dz  , 


'th 


2th  ° 


by  a  change  of  variables 

w  =  z2  dw  =  2zdz 


W  =  t L 
"th  th 


then 


PFA  =  fz  ~  exp[-w/2a2]dw  =  exp[-Z2h/2o2]  .  (2.4) 

Zth  2a 

This  equation  can  be  used  to  determine  a  threshold  given  a  desired 
probability  of  false  alarm,  i.e., 

(2.5) 


_2A 


Zth  =  [-ln(PFA)2otr 
2.1.2  Square  Law  Detector 

A  square  law  detector  produces  an  output  which  is  proportional  to 
the  square  of  the  video  envelope  and  is  given  by 

(2.6) 


2  2 

y  =  xl  +  xQ  , 


where  Xj  is  the  in-phase  video,  Xq  is  the  quadrature  video,  and  y  is 
the  detector  output. 

If  Xj  and  Xq  are  independent  zero  mean  Gaussian  random  variables 

2 

and  homogeneous,  i.e.,  they  have  the  same  variance,  a  ,  the  pdf  of 
y  is 


P(y)  =  -V  exp[-y/2a2]  . 

2<j 

The  pdf  is  the  well  known  exponential  distribution. 


(2.7) 


Assuming  a  fixed  threshold,  Y^,  the  probability  of  false  alarm  is 

/co 

p(y)dy  =  exp[-Yth/2o2J  .  (2.8) 

Yth 

This  equation  can  be  used  to  determine  a  threshold  given  a  desired 
probability  of  false  alarm,  i.e., 

Yth  =  [-1 n ( PFA) 2a2]  .  (2.9) 

2.2  Target  Plus  Noise 

This  section  gives  the  pdf  and  probability  of  detection  expressions 
for  the  single  pulse  amplitude  detected  target  plus  noise  cases.  The 
target  will  be  either  a  nonfluctuating  or  steady  target  or  a  Swerling  I 
target. 

At  this  point,  it  is  desirable  to  discuss  the  definition  of  inter¬ 
mediate  frequency  (IF)  signal -to-noise  ratio,  x,  commonly  found  in  radar 
literature.  The  basic  writings  of  Marcum  [1],  Swerling  [2],  and  Rice  [3] 
used  the  following: 

2 

Average  Signal  Power  at  IF  _  JP _  in\ 

x  Average  Noise  Power  at  IF  ^2  '  ’ 

where  the  received  target  is  Pcos(2irft  +  6)  whose  IF  average  power  is 
2 

P  /2. 

For  simplicity,  it  can  be  assumed  that  any  quadrature  channel 
processing,  such  as  a  clutter  rejection  filter,  will  not  affect  the 
signal -to-noise  ratio.  Hence  the  IF  signal -to-noise  ratio  and  the 
detector  input  signal -to-noise  ratio  are  the  same. 

2.2.1  Steady  Target 

A  steady  target  [1]  is  defined  as  a  target  where  the  signal -to- 


noise  ratio  for  one  pulse  describes  the  signal -to-noise  ratio  of  any 
pulse  of  a  train  under  consideration.  Hence,  the  pdf  is 
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p(x)  =  P  ‘'/2a  ,  (2.11) 

where  P  is  the  IF  signal  peak. 

2. 2. 1.1  Linear  Detector 

The  pdf  of  a  single  steady  target  plus  noise  variate,  z,  after 
linear  detection  is  derived  in  Appendix  A: 


p(z)  =  \  exp 


z2  +  P21  T  /zP\ 

L'^Tn  “la2/’ 


(2.12) 


2 

where  P  is  the  peak  signal  voltage  prior  to  detection,  a  is  again  the 

same  variance,  and  I  is  the  modified  Bessel  function  of  the  first  kind 

o 

of  zero  order  [4]. 

The  probability  of  detection  is  given  by 

oo  oo  "**  2  2  "1 

PD  =  /  p(z)dz  -  /  4  exp  -  z  I  (^)dz  ■  (2-13) 

Zth  Zth°  L  2°  J  W 


This  equation  is  of  the  form  of  a  Q-function  [4] 


oo  r  2  2~ 

Q(b,c)  =  ^  a  exp  -  5-+-L.  ij 

and  thus  Equation  (2.13)  can  be  written  as 
PD  =  Q(b,c)  , 


IQ(ab)da 


(2.14) 


(2.15) 


where  b  =  P/a  and  c  =  Z^./o. 

th 

2 . 2 . 1 . 2  Square  Law  Detector 

The  pdf  of  a  single  steady  target  plus  noise  variate,  y,  of  a 


square  law  amplitude  detector  is 


p(y,?exp["^]  '°(^) 


(2.16) 


which  is  obtained  from  Equation  (2.12)  by  a  change  of  variables. 
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The  probability  of  detection  for  the  square  law  detector  is  the 
same  as  for  the  linear  detector  given  in  Equation  (2.15),  i.e., 

PD  =  Q(b,c)  ,  (2.17) 

where  b  =  P/a  and  c 

2.2.2  Swerling  I  Target 

A  Swerling  I  target  [2]  is  defined  as  samples  which  are  correlated 
within  a  pulse  group  but  are  independent  on  a  scan-to-scan  basis 
(slowly  fading).  This  case  is  applicable  to  many  radar  targets  since 
they  tend  not  to  be  independent  from  pulse  to  pulse,  but  independent 
from  scan  to  scan  due  to  target  position  change. 

The  pdf  for  a  single  sample  signal -to-noise  ratio,  x,  is 

w(x,x)  =  ~  expT-  *  (2.18) 

—  ~2  2 

where  x  =  p  /2a  is  the  average  signal -to-noise  ratio. 


2. 2. 2.1  Linear  Detector 

The  mathematical  analysis  of  a  Swerling  I  target  plus  noise  and 
a  linear  detector  is  difficult  and  no  analysis  was  found  in  the 
literature. 

2 . 2 . 2 . 2  Square  Law  Detector 

The  pdf  for  a  single  Swerling  I  target  plus  noise  square  law 
detector  output  variate,  y,  is  derived  in  Appendix  A  and  is  given  as 


p(y)  = 


i 

2 a2  (1  +x) 


X _ ~ 

(1  +50. 


u(y)  • 


(2.19) 


The  probability  of  detection  is  given  by 


Substituting  Equation  (2.8)  yields 

1 

PD  =  PFA  1  +  *  .  (2.21) 

2.3  Concl usi ons 

A  fixed  threshold  decision  element  is  normally  used  to  specify 
radar  system  performance.  Due  to  the  complex  equations  obtained  when 
a  linear  detector  and/or  a  steady  target  is  used,  the  performance  will 
normally  be  based  on  a  Swerling  I  target  model  and  a  square  law 
detector.  This  assumption  does  not  cause  any  significant  problems. 

The  Swerling  I  target  model  is  a  realistic  model  for  many  radar  targets 
and  the  square  law  and  linear  detector  have,  as  shown  by  Marcum  [1], 
essentially  the  same  detection  performance  for  a  single  pulse. 

In  an  actual  radar  system,  the  use  of  a  fixed  threshold  would 

.  p> 

require  having  a  priori  knowledge  of  the  thermal  noise  variance,  a  , 
to  maintain  a  desired  probability  of  false  alarm.  For  example,  if 
PFA  =  10'^,  then  from  Equation  (2.8) 

Yth  =  -ln(PFA)2o2  =  27.63a2  .  (2.22) 


2 

Hence  is  a  function  of  the  input  noise  variance,  a  .  As  shown  in 
Chapter  I,  the  probability  of  false  alarm  is  strongly  affected  by  a 
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Even  if  exact  knowledge  of  the  thermal  noise  were  available,  the 
total  system  interference  variance  can  change  due  to  residual  clutter 
not  cancelled  by  the  prewhitening  filter  or  jammers.  Therefore,  an 
adaptive  technique  for  determining  the  threshold  is  required.  These 
techniques  are  referred  to  as  constant  false  alarm  rate  (CFAR) 
processors  or  adaptive  detection  processors. 


CHAPTER  III.  ADAPTIVE  THRESHOLD  ANALYSIS 


3.0  Introduction 

This  chapter  reviews  the  theoretical  analysis  of  three  commonly 
found  CFAR  processors:  cell  averaging,  "greatest-of and  log. 
Basically,  these  processors  sample  the  background  interference  in  the 
time  domain  around  a  cell,  i.e.,  a  range  cell  of  interest,  and  then 
utilize  the  samples  to  estimate  the  unknown  statistical  parameters  of 
the  interference.  This  estimate  is  used  to  determine  a  threshold  for 
the  cell  of  interest. 

The  estimated  threshold's  probability  density  functions  are  given, 
and  equations  for  the  probability  of  false  alarm  and  probability  of 
detection  are  derived  for  the  cell  averaging  and  "greatest-of"  CFAR 
techniques.  The  analysis  assumes  a  square  law  detector  and  a  Swerling 
I  target  for  reasons  stated  in  Chapter  II. 

Only  a  limited  analysis  of  the  log  CFAR  is  presented  due  to  a  lack 
of  available  analytical  results.  An  equivalence  to  the  cell  averaging 
technique  is  discussed. 

As  in  the  fixed  threshold  analysis,  the  CFAR  processor  analysis 
will  be  based  on  white  Gaussian  noise  interference  which  is  a  result 
of  the  prewhitening  or  clutter  rejection  filter. 

3.1  Cell  Averaging  CFAR  Analysis 

This  procedure  (Figure  3)  forms  the  threshold  Yth  by  scaling  the 
average  value  of  N  square  law  detected  reference  cell  outputs  of  the 
quadrature  channels,  I  and  Q,  i.e.. 
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y  N  v  N  9  „  2N 

Yth  =  T  yn  =  T  (XI  +  xC))n  =  T  nS,  ^xIQ^n 

(3.1) 

where  the  last  summation  results  since  it  is  equivalent  to  summing  2N 

statistically  independent,  squared,  zero  mean  Gaussian  random  variables 

Xjq.  It  is  assumed  that  the  referenced  cells  are  homogeneous,  i.e., 

2 

each  (xjq)n  has  the  same  variance,  a  .  Consequently  the  distribution 
2 

for  NY^/Ka  will  have  a  chi-square  pdf  with  2N  degrees  of  freedom. 

The  pdf  for  is  obtained  by  changing  variables  on  the  chi-square 
pdf  [5].  Thus 


Figure  3.  Block  Diagram  of  a  Conventional 
Cell  Averaging  CFAR  Processor 


This  equation  can  be  used  with  the  fixed  threshold  PFA  of  Equation 
(2.8)  to  obtain  the  expected  PFA  when  the  ceil  averaging  CFAR  controls 
the  threshold,  i.e.. 
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f  ' 

i**» 


h 

-  * 


PFACft  =  EiPFA)  -jf  exp[-Vt|)/2oJJp(Vth)dYth 

■  (if  wAn  ^  r  ex(,[-vt^/2a2](^)N  exp[-1 


2o  0 


NYth/2Ka  IciV 


th 


Letting  a  =  Yt^/(2o^)  gives 


PFACA  =  (1)N  liTTTr  /  aN_1  exp[-a(N/K+l)]da  .  (3.3) 


Letting  a  =  b/(N/K+l )  yields 

vN 


pfaca  ■  (l)  unrr  I 


°gh 


■  N-l  Wk 

exp[-b] 


N-l 


N  +  1 

K 


=  (I)'  WOT  T^N  f0  bN'’  «Pt-bJ«b 


pfaca  ■  (’  +  Ij"  ■ 


This  allows  the  CFAR  threshold  constant,  K,  to  be  determined  from 
the  desired  average  probability  of  false  alarm,  PFA^,  i.e., 


(3.4) 


It  is  easily  seen  that  the  average  probability  of  false  alarm  is 
not  dependent  on  the  noise  variance.  Hence,  the  Gaussian  noise  level 
does  not  have  to  be  known  to  maintain  CFAR.  Nitzberg  [6]  called  this 
an  unknown  level  CFAR,  but  it  is  commonly  known  as  a  range  cell 
averaging  CFAR. 

The  expected  value  of  the  probability  of  detection  PD^  for  the 
Swerling  I  target  can  be  determined  by  the  same  procedure,  i.e.. 
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P0CA  =  E{PO}  -  J  PO  P(Yth)dYth 

—OO 


£exp[-Yth/202(l«)]p(vth)dYth 


(3.5) 


Substituting  Equation  (3.3)  yields 


1  4-  x 

*  + 


-1/N 


N 


where  x  is  the  average  IF  signal-to-noise  ratio. 

This  result  can  be  used  to  plot  PDCA  versus  "x  with  PFA^  and  N  as 
parameters.  A  typical  curve  is  shown  in  Figure  4. 
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Since  the  fixed  threshold  performance  curves  are  extensively 
tabulated  [7],  a  general  cell  averaging  CFAR  signal -to-noise  loss  curve 
is  desirable.  The  loss  is  given  by  the  equation  from  Moore  [8], 


L  = 


-10  log 


r  -r/n 

-  1 

,  1  +  XCA 

9  x  +  10R/N 

XCA  IU  J 

XCA 

(3.6) 


-R 


where  R  corresponds  to  the  exponential  in  the  PFA^,  i.e.,  PFA^  -  10  , 


CA 


N  is  the  number  of  reference  cells,  and  x^  is  cell  averaging  signal-to- 
noise  ratio  necessary  to  give  the  same  PD  at  x  for  a  fixed  threshold 
detector.  This  loss  is  shown  in  Figure  5. 


Figure  5.  Cell  Averaging  CFAR  Loss 


3.2  "Greatest-Of"  CFAR  Analysis 

In  this  method  (Figure  6),  the  reference  cells  are  divided  into 
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two  subsets  of  size  N/2.  The  cell  averaging  method  is  used  to  deter¬ 
mine  a  threshold  for  each  of  the  reference  cell  subsets.  One  subset 
is  located  before  the  reference  cell  of  interest  and  the  other  after 
the  reference  cell  of  interest.  The  "greatest-of"  CFAR  threshold  is 
obtained  by  selecting  the  largest  value  from  the  two  subset 
thresholds,  i.e., 

Kr  N/2 

Y1  =  N/2  yn 
1C  -N/2 

Y2  =  N/2  njC  1  yn 

Yg  »  MAX  [Yr  Y2]  .  (3.7) 

The  G  subscripts  for  "greatest-of"  are  used  so  that  there  is  a 
distinction  from  the  cell  averaging  processor. 


Figure  6.  Block  Diagram  of  a  "Greatest-of"  CFAR  Processor 
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Since  Y-j  and  Y£  are  determined  by  the  cell  averaging  method,  then 
2  2 

NY-|/(2KgO  )  and  NY2/(2Kgcr  )  are  chi-square  distributions  with  N  degrees 
of  freedom.  Consequently,  the  pdf  for  Y-j  (or  Y2)  can  be  obtained  from 
Equation  (3.2). 

By  replacing  Y  with  Y-|  (or  Y2)  and  N  with  N/2,  i.e.,  for  Y-j 


plv,)  =  ^ 


N/2  -  1 


exp 


N  V 

■a?  57 


u(Y1). 

(3.8) 


Papoulis  [9]  gives  an  expression  for  finding  a  pdf  of  the  maximum  of 
two  random  variables,  cf..  Equation  (7-15),  p.  193, 

PG(VG)  =  2F(Y)p(Y)|¥=Yg 

•  2F¥(Yg)fv(Yg)  (3.9) 

where  the  cumulative  distribution  function  for  Y1  and  Y2  is  represented 
by  F(-). 

The  average  probability  of  false  alarm  for  the  "greatest-of"  CFAR 
is  derived  in  Appendix  A  and  is  given  by 


2  PFA_ 


pfag  = 


(n  +  h')' 

(I  -  1)1(1  +  PFAp'2/Ny/2  0=0  n !  (l  +  PFAp2/N) 


N/2  -  1 
1  N/2  E. 


(3.10) 


where  PFAp  is  called  the  prototype  section  PFA  and  is  equal  to 

-N/2 


PFA„ 


(3.11) 


Note  that  this  is  an  expression  for  a  cell  averaging  CFAR  which  uses 
N/2  reference  cells. 


It  is  possible  to  solve  for  Kg  in  terms  of  PFAp  as  in  Equation 
(3.4)  with  N  replaced  by  N/2,  but  Kg  is  not  easily  related  to  PFAg 
However,  the  results  of  Equation  (3.10)  can  be  plotted  as  shown  in 
Figure  7,  then  used  to  obtain  the  threshold  constant.  For  example,  if 
it  is  desired  to  establish  PFAg  =  10"^  with  N  =  32,  then  from  Figure  7, 

_  _C 

PFAp  =  1.75  x  10  .  Consequently,  Kg  is  calculated  to  be  15.73  and 

would  be  used  in  the  Y-|  and  determinations  in  order  to  establish 
PfAg  =  10'6. 


Figure  7.  GO  CFAR  False  Alarm  Characteristics 

Similar  results  derived  in  Appendix  A  hold  for  the  probability  of 
detection 
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PD, 


-  /;  P°  PG<VdVC 


2  PDr 


(I  -  l)  !  (l  t  P0p-2/N) 


N/2-1  (n  +  ^  -  l)l 

N/2  L 


n=0 


n!( 


1  +  PDp_2/N) 


n 


where  the  prototype  section  PDp  is 

-N/2 


2Kr 


x  +  PFA, 


-2/N 


-N/2 


PDp  =  [1  + 


NO  +  x) 


1  +  x 


(3.12) 
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This  represents  the  performance  of  a  cell  averaging  CFAR  with  N/2 
reference  cells.  Once  PFAp  has  been  found  (as  from  Figure  7),  then  PDp 
can  be  calculated  from  Equation  (3.13)  and  PDg  from  Equation  (3.12). 

It  would  be  highly  desirable  to  determine  the  signal -to-noise  ratio 
loss  for  the  "greatest-of "  CFAR  as  compared  to  the  ideal  fixed  threshold. 
Unfortunately,  the  complexity  of  Equation  (3.12)  prevents  such  an 
analysis. 

Analysis  of  the  "greatest-of"  CFAR  has  been  performed  [8,  10,  11, 
13].  One  advantage  of  the  "greatest-of"  CFAR  is  discussed  in  References 
8,  10,  and  11,  that  is,  the  improved  regulation  of  false  alarms  obtained 
for  range  extended  clutter  when  compared  to  a  cell  averaging  CFAR. 

Range  extended  clutter,  discussed  further  in  Chapters  IV  and  V,  is  the 
weather  or  chaff  clutter  not  rejected  by  the  clutter  filter  and  occupy¬ 
ing  some  of  the  CFAR  reference  cells. 

3.3  Log  CFAR  Analysis 

This  system  (Figure  8)  forms  an  estimate  for  the  threshold  as 


VT  = 


K  n  y.1/N=K 
j,l  J 


N 

n  y, 

LH  JJ 


1/N 


(3.14) 


An  equivalent  method.  Figure  9,  for  processing  is  to  use  a  log 
detector  at  the  input  such  that 

N  , 

log  VT  =  £  n-  log  y.  +  log  K  (3.15 

1  j=l  N  J 

is  formed  and  log  yQ  is  compared  to  this  threshold. 

The  expected  value  of  the  estimate  is  determined  to  be 

eivti  .  4[yn)N ■ 

The  gamma  function  will  become  approximately  equal  to  1  for  large 
N  since  T(l)=l.  Thus  the  expected  value  of  the  threshold  will  become 


ih  i) 


(3.16 


LIM  E{VT}  =  LIM  2K a2 

oo  N  -*■  CO  j 


(3.17 
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which  is  equal  to  2KaQ  for  homogeneous  noise.  Therefore  a  reasonable 
estimate  can  be  formed  by  using  the  log  CFAR  algorithm. 

Whereas  a  detailed  mathematical  analysis  has  not  been  performed, 
Hansen  and  Ward  [12]  have  performed  a  Monte  Carlo  analysis  of  the  log 
CFAR.  Nitzberg  [6],  concerning  a  similar  algorithm  called  the 
geometric-mean  CFAR,  has  determined  the  probability  of  detection  when 
an  assumption  is  made  about  the  noise  distribution  in  the  auxiliary 
cells,  viz.,  the  geometric-mean  assumption. 

In  comparing  the  log  CFAR  and  the  cell  averaging  CFAR,  Hansen  and 
Ward  [12]  have  proposed  an  empirically  determined  formula  for  the 
relationship  between  the  number  of  reference  samples  required  by  the 
two  detectors  in  order  for  their  CFAR  losses  to  be  identical : 

Nlog  =  1-65  NCA  -  °-65  •  (3J8 

The  main  advantage  of  the  log  CFAR  is  the  increased  dynamic  range 


available  due  to  the  log  detector. 


3.4  Summary 

The  cell  averaging  CFAR  has  been  utilized  extensively  in  radar 
signal  processors  due  to  its  capability  in  homogeneous  noise  and  its 
well  understood  and  analyzed  performance. 

The  log  CFAR  is  simply  a  cell  averaging  CFAR  following  a  log 
detector  which  provides  performance  equivalent  to  the  cell  averaging 
CFAR  if  the  number  of  cells  is  sufficient.  The  log  CFAR  has  been  used 
extensively  due  to  its  dynamic  range  capability. 

The  "greatest-of "  CFAR  has  not  been  used  extensively,  due  partly 
to  the  original  belief  that  it  had  approximately  a  1  dB  loss  over  the 
cell  averaging  CFAR,  e.g.,  Hansen  [10].  Recent  work  by  Moore  [8], 
Moore  and  Lawrence  [11]  and  Hansen  and  Sawyer  [13]  has  shown  only  a 
0.2  dB  difference  in  the  two  processors.  Hence,  the  "greatest-of" 
CFAR,  whose  main  advantage  is  the  improved  false  alarm  regulation  in 
extended  clutter  [8,  11]  should  have  increased  utilization  in  radar 
processors. 


CHAPTER  IV.  DESCRIPTION  OF  SIMULATION 


4.0  Introduction 

The  cell  averaging  CFAR  and  the  "greatest-of"  CFAR  are  two  commonly 
used  techniques.  Analysis  of  the  cell  averaging  CFAR  has  been  exten¬ 
sively  performed  [10,  14-16 J ,  but  only  limited  analysis  of  the  "greatest- 
of"  CFAR  has  been  performed  [8,  11,  13]. 

The  main  thrust  of  this  study  is  to  determine  the  performance  of 
the  two  CFAR  processors  by  development  of  a  simulation  and  utilization 
of  Monte  Carlo  techniques.  The  performance  results  obtained  are  used  to 
compare  the  two  techniques.  This  chapter  gives  a  description  of  the 
simulation. 

4.1  Simulation  Description 

A  block  diagram  of  the  simulation  is  shown  in  Figure  10. 


Figure  10.  Block  Diagram  of  CFAR  Simulation 


First,  synthetic  video  composed  of  a  combination  of  target,  noise 
and  clutter  is  generated  for  a  quadrature  channel  processor.  The 
amplitude  is  extracted  by  an  exact  square  law  or  linear  detector.  The 
detected  output  is  compared  against  a  threshold  determined  by  either  a 
cell  averaging  CFAR  or  a  "greatest-of"  CFAR  using  other  detected 
outputs . 
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If  the  detector  output  exceeds  the  threshold,  a  target  detection 
is  reported.  If  the  output  is  interference  only,  this  is  a  false 
alarm;  if  the  output  contains  signal,  this  is  a  detection.  A  detection 
and  false  alarm  count  are  maintained  for  both  processors.  Finally, 
after  a  number  of  Monte  Carlo  trials  the  detection  and  false  alarm 
counts  are  used  to  calculate  a  probability  of  detection  and  a  proba¬ 
bility  of  false  alarm  for  each  processor. 

4.2  Synthetic  Video 

This  section  discusses  the  target  models,  noise,  and  clutter  used 
in  the  simulation. 

4.2.1  Target  Models 

Two  target  models  were  used  in  the  simulation:  a  steady  or  non¬ 
fluctuating  target  [1]  and  a  Swerling  I  target  [2]. 

The  steady  target  is  defined  as  a  target  where  the  signal -to-noise 
ratio  for  one  pulse  describes  the  signal -to-noise  ratio  of  any  pulse  of 
a  train  under  consideration.  A  steady  target  is  modeled  in  the  I  and 
Q  channels  by 

Sj  =  Pco$(0) 

Sq  =  Psin(G)  (4.1 ) 

where  P  is  the  IF  peak  signal  voltage  and  9  is  a  uniformly  distributed 
random  phase  angle. 

A  Swerling  I  target  is  defined  as  samples  which  are  correlated 
within  a  pulse  train  but  are  independent  on  a  scan-to-scan  basis 
(slowly  fading).  This  case  is  applicable  to  many  radar  targets,  since 
they  tend  not  to  be  independent  from  pulse  to  pulse,  but  due  to  target 
position  change,  independent  from  scan  to  scan.  The  probability 
density  function  for  one  sample,  x,  is 
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where  x  =  P  /2a  is  the  average  signal -to-noise  ratio.  Since  the  power 
distribution  of  a  Swerling  I  target  is  the  well  known  exponential,  then 
the  amplitude  distribution  is  Rayleigh  and  the  Swerling  I  target  models 
in  the  I  and  Q  channels  are  given  by 

S j  =  P  ^2Tnu-|  cos(2ttu2) 

Sg  =  P  ^2TntTj  sin(2ru2)  (4.3) 

where  u-j  and  u2  are  independent  uniformly  distributed  variates  from 
0  to  1 . 

4.2.2  Noise  Model 

The  system  noise  will  be  zero  mean  Gaussian  noise  whose  pdf  is 
given  by 

p(v)  =  — - —  exp[-v^/2a^]  (4.4) 

V27 ra2 

2 

where  a  is  the  variance. 

There  is  a  procedure  for  generating  uncorrelated  Gaussian  samples 
called  the  direct  method  [17].  In  this  procedure,  pairs  of  independent 
samples  (u-j ,  u2)  are  drawn  from  a  uniform  distribution  (0  to  1),  then 
transformed  as 

v-j  =  -^-21  nu^  cos(2ttu2) 

v2  =  y-21nu2  sin(2iru2)  (4.5) 

where  v^ ,  v2  are  the  uncorrelated  samples  of  the  Gaussian  distribution. 
From  Equation  (4.5),  and  u2  may  be  expressed  as  functions  of  v^ 
and  v2 
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(4.13) 


The  variance  of  v-|  and  is 


Hence,  the  Gaussian  noise  is  modeled  in  the  I  and  Q  channels  as 
Ni  =  ovi 

Nq  =  av2  (4.14) 

where  v-|  and  are  defined  in  Equation  (4.5)  and  a  is  the  standard 
deviation  in  each  channel  and  at  IF. 

4.2.3  Clutter  Model s 

Two  clutter  models  were  included  in  the  simulation:  nonhomogeneous 
interference  and  Weibull  [18]  distributed  clutter. 

The  nonhomogeneous  interference  is  clutter  where  the  power  density 
varies  as  a  function  of  range,  i.e.,  chaff  or  weather  clutter  which  is 
distributed  in  range.  The  clutter  power  appears  as  a  step  function 
with  a  clutter  edge  [14]  as  shown  in  Figure  11. 


Figure  11.  Clutter  Edge  Model 


cl 


The  clutter  is  assumed  to  be  Gaussian  in  each  range  cell  and  the 
utter  powers  in  the  N  CFAR  reference  cells  are  related  by  a  ratio 
such  that 


T_  = 


n  =  1  ,  2,  ...  N 


(4.15) 


■  **~*X- 
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2  2 
where  o  is  the  clutter  power  in  the  cell  of  interest  and  an  is  the 

clutter  power  in  the  reference  cell. 

In  Reference  18,  Boothe  has  shown  that  the  spatial  distribution  of 

the  ground  clutter  backscatter  coefficient,  o°,  for  various  types  of 

terrain  fit  quite  well  with  a  Weibull  pdf.  The  Weibull  pdf  is  given  by 


where  b  =  1/A  (A  =  Weibull  slope  parameter)  and 


(4.16) 


a  = 


1  n2 


(4.17) 


where  =  median  value  of  Weibull  pdf.  Typical  values  of  the  clutter 
slope  parameter  (A)  and  median  backscatter  coefficient  (o°)  are  given 
in  Reference  18. 


A  single  Weibull  sample,  a ,  can  be  generated  by 
o 


o 

o  = 


(ln2  Y 


[-ln(u)] 


(4.18) 


where  u  is  a  uniformly  distributed  random  variate.  Due  to  the  quadra¬ 
ture  channel  processing,  two  independent  Weibull  samples,  oj  and  Oq, 
must  be  generated. 

Hence,  the  Weibull  pdf  is  given  in  the  I  and  Q  channels  as 


(ln2) 


A  [-ln(u1 )]  cos(2ttu2) 


o°  -^hr  ['ln(u,)]A  sin(2™?)  (4.19) 

w  (ln2)ft  1  c 

where  u-|  and  u2  are  independent  samples  drawn  from  a  uniform  distribu¬ 
tion  (0,1 ) . 


4.3  Detector  Laws 

The  square  law  detector  output  is  given  by 

y  =  I2  +  Q2  h  (4.20) 

where  I  and  Q  are  the  video  in  the  in-phase  and  quadrature  channels, 
i.e.,  signal  plus  interference,  respecti vely. 

The  linear  detector  output  is  given  by 


where  I  and  Q  are  as  above. 

4.4  CFAR  Processors 

Two  constant  false  alarm  rate  processors  are  modeled  in  the 
simulation. 


where  the  yn's  are  the  detector  outputs,  n  is  the  reference  cell  index 
and  K  is  the  scaling  constant.  The  actual  model  is  implemented  as 
shown  in  Figure  3,  that  is, 

u  r-N/2  N/2 

Yth'T  E.  <4-23 

n=  - 1  n= I 

where  yQ,  the  cell  of  interest,  is  not  included  in  the  threshold 


determination. 


4.4.2  "Greatest-Of"  CFAR 


The  "greatest-of"  CFAR  will  form  a  threshold  Yg  by  using  the  cell 
averaging  CFAR  processor  on  two  sets  of  N/2  detected  reference  cell 
outputs  and  will  select  the  largest  value  obtained.  The  "greatest-of" 
processor  is  simulated  as 


Kr  -N/2 

Y1  =  N72  nP_}  yn 


Kr  N/2 

Y2  =  N72  ^  yn 


(4.24) 


and 

Yg  =  MAX  [Yr  Y2] 

where  Kg  is  the  "greatest-of"  scaling  constant.  Again  y  ,  the  cell  of 
interest,  is  not  included  in  the  threshold  determination. 

4.5  PFA  and  PD  Determinations 

The  probabilities  of  false  alarm  are  determined  from  detector 
outputs  which  contain  noise  and/or  interference  only.  A  threshold  for 
the  cell  of  interest  is  calculated  b,  uhe  cell  averaging  CFAR  processor 
using  other  detector  outputs.  The  magnitude  of  the  cell  of  interest 
is  compared  to  the  threshold  and,  if  it  is  larger,  a  false  alarm  is 
reported.  In  the  simulation,  a  false  alarm  counter  (FAC)  is  initial¬ 
ized  to  zero  at  the  beginning  of  a  Monte  Carlo  sequence,  then  FAC  is 
incremented  by  one  for  each  false  alarm  reported.  Finally,  the 
average  probability  of  false  alarm  is  determined  as 


FAC 

NMON 


(4.25) 


where  NMON  is  the  number  of  Monte  Carlo  trials.  The  average  probability 
of  false  alarm,  PFAg,  for  the  "greatest-of"  CFAR  is  determined  by  the 
same  procedure. 
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The  probabilities  of  detection  are  determined  from  detector  outputs 
which  contain  targets.  The  same  procedure  is  used  as  in  the  PFA  deter¬ 
mination;  however,  a  target  detection  counter  (TDC)  is  incremented  for 
each  threshold  that  is  exceeded  by  the  magnitude  of  the  cell  of  interest 
in  which  a  target  resides.  Then  the  average  probability  of  detection 
is  determined  as 


PD 


TDC 

CA  ~  NMON  ' 


(4.26) 


The  average  probability  of  detection,  PD^ ,  for  the  "greatest-of" 
is  determined  in  a  similar  manner. 

4.6  Summary 

A  simulation  has  been  developed  which  can  be  used  to  determine 
the  cell  averaging  CFAR  and  "greatest-of"  CFAR  performance  for  different 
environmental  conditions,  targets  and  detectors.  The  probability  of 
false  alarm  and  the  probability  of  detection  results  obtained  can  be 
used  to  verify  the  theoretical  performance  equations  and  to  compare  the 
relative  performance  of  the  two  processors. 
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CHAPTER  V.  RESULTS 


5.0  Introduction 

The  simulation  described  in  Chapter  IV  was  developed  to  compare 
the  performance  of  the  cell  averaging  and  "greatest-of "  CFAR  processors. 
The  utilization  of  a  simulation  allows  determination  of  the  processors' 
performance  for  the  different  targets,  detectors  and  clutter  environ¬ 
ments  simulated.  The  probabilities  of  false  alarm  and  probabilities  of 
detection  are  the  basis  for  comparing  the  two  CFAR  processors.  The 

desired  probability  of  false  alarm  in  radars  is  normally  quite  small, 

-3  -9 

i.e.,  10  to  10  .  Thus,  it  is  difficult  to  verify  the  probability 

of  false  alarm  using  a  computer  simulation  due  to  the  amount  of 

computer  time  required  to  complete  a  sufficient  number  of  Monte  Carlo 
5 

passes,  i.e.,  10  or  more.  This  difficulty  was  overcome  by  programming 
the  simulation  on  an  array  processor.  The  array  processor  is  a  high 
speed  arithmetic  unit  designed  for  scientific  applications.  A  brief 
discussion  of  the  array  processor  is  given  in  Appendix  D. 

5.1  Probability  of  False  Alarm  Results 

To  compare  the  two  CFAR  techniques  it  is  necessary  to  design  them 
to  maintain  the  sane  average  probability  of  false  alarm  in  homogeneous 
noise,  i.e.,  PFAq  =  PFA^. 

For  the  cell  averaging  CFAR  it  is  only  required  to  specify  the 
desired  average  probability  of  false  alarm,  PFAC^,  and  the  number  of 
reference  cells  N,  and  by  using  Equation  (3.4)  to  determine  the 
threshold  constant  K. 
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For  the  "greatest-of"  CFAR  the  design  procedure  is  somewhat 


complicated.  The  threshold  constant  Kg  is  determined  by  a  computer 
program  which  iterates  PFAp  in  Equation  (3.10)  until  the  desired  value 
for  the  PFAg  is  obtained.  Then  this  value  of  PFAp  and  the  number  of 
reference  cells  N  is  used  in  Equation  (3.11)  to  determine  Kg. 

Hence,  theoretically  PFAgA  =  PFAg  Tor  the  same  number  of  reference 
cells  in  homogeneous  noise. 

-3  -4  -5 

The  design  probabilities  of  false  alarm  of  10  ,10  ,  and  10 

were  chosen  because  they  are  commonly  found  values  and  because  the 
Gaussian  random  number  generator  lacks  distribution  tails  necessary  for 

_5 

a  false  alarm  rate  <10  .  The  CFAR  window  widths  N  were  chosen  to  be 

8,  16,  and  32,  since  digital  hardware  is  normally  implemented  in  powers 
of  two.  The  probabilities  of  false  alarm  for  the  cell  averaging  and 
"greatest-of"  CFAR  obtained  by  the  Monte  Carlo  simulation  are  given  in 
Table  1. 


Table  1.  CFAR  Processor  Probabilities  of  False  Alarm 


PFA0 

N 

PFAca 

PFAp 

PFAg 

8 

0.113-2 

0.877-2 

0.106-2 

10"3 

16 

0.110-2 

0.507-2 

0.107-2 

32 

0.104-2 

0.318-2 

0.103-2 

_A 

8 

0.121-3 

0.223-2 

0.123-3 

10"4 

16 

0.104-3 

0.104-2 

0.110-3 

32 

0.110-3 

0.532-3 

0.124-3 

c 

8 

0.106-4 

0.599-3 

0.770-5 

10"5 

16 

0.134-4 

0.230-3 

0.144-4 

32 

0.115-4 

0.940-4 

0.115-4 

The  number  of  Monte  Carlo  runs  used  to  determine  the  results  in 
5  6  7 

Table  1  were  10  ,  1C  ,  and  10  for  the  probabilities  of  false  alarm 


i 
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10  ,  10  ,  and  10  ,  respectively.  The  number  of  Monte  Carlo  runs 
required  to  give  a  priori  probabilities  PFA  and  PD  for  a  specified 
range  of  the  estimated  parameters  is  calculated  in  Appendix  C. 

5.2  Probability  of  Detection  Results 


The  performance  curves  in  Figures  12  through  23  were  determined  by 
the  Monte  Carlo  simulation.  The  curves  are  plotted  as  probability  of 
detection  versus  input  signal -to-noise  ratio.  The  cell  averaging  CFAR 
performance  curves  are  given  in  Figures  12  through  14  for  a  steady 
target  and  in  Figures  15  through  17  for  a  Swerling  I  target.  A  square 
law  detector  is  used.  The  "greatest-of "  CFAR  performance  curves  are 
given  in  Figures  18  through  20  for  a  steady  target  and  in  Figures  21 
through  23  for  a  Swerling  I  target.  Again,  a  square  law  detector  is 
used. 
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The  design  false  alarm  probabilities  of  10  ,10  ,  and  10  are 


shown  on  the  plots  while  the  actual  average  false  alarm  probabilities 
are  given  in  Table  1 . 

The  theoretical  probability  of  detection  equations  have  been 
determined  for  the  ceil  averaging  CFAR  with  a  steady  target  and  a 
Swerling  I  target  and  for  the  "greatest-of"  CFAR  with  a  Swerling  I 
target  only.  These  equations  are  shown  below. 

The  theoretical  PD  for  a  square  law  detected  steady  target  and 
cell  averaging  CFAR  is  derived  in  Reference  14  and  is  given  as 

2  N-l  _ r  -2W_2 

r..cn.  f 

g~  +  i  m=o  ""  'g' 

where  g  =  Ky2/N  is  the  input  signal -to-noise  ratio,  e  =  a  g  /[ 2(g  +2)] 
and  L  ( e )  are  Laguerre  polynomials  with  the  properties 


PD  =  1  - 


i!_  V  sat:  /-/iainu  (2  y  L  (e).  (5.D 

+  1  mo  n!  V  *  2/  m 


L  (e)  =  1,  L,(e)  =  1  +  e,  L  Me)  =  (e  +  2m+  1  )L  (e)  -  nrL  ,  (e). 
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Figure  12.  Cell  Averaging  CFAR  Performance  Curves  (PFA, 
Steady  Target,  Square  Law  Detector)  1 


SNR  IdB) 

Averaging  CFAR  Performance  Curves 
eady  Target,  Square  Law  Detector) 


Cell  Averaging  CFAR  Performance  Curves 
Swerling  I  Target,  Square  Law  Detector) 
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Figure  18.  "Greatest-Of "  Performance  Curves  (PFA, 
Steady  Target,  Square  I  *w  Detector) 
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Figure  19.  "Greatest-Of"  Performance  Curves  ( 
Steady  Target,  Square  Law  Detector) 
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Figure  20.  "Greatest-Of"  Performance  Curves  (PFAD  =  10  , 

Steady  Target,  Square  Law  Detector) 
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The  theoretical  PD  for  a  square  law  detected  Swerling  I  target 
using  cell  averaging  CFAR  is  given  in  Equation  (3.5)  and  repeated  here. 


The  probability  of  detection  curves  for  an  ideal  threshold,  i.e., 
a  fixed  threshold  system  where  the  noise  power  is  known,  are  given  in 
Meyer  and  Mayer  [7].  Ideal  threshold  curves  are  plotted  in  Figures  12, 
15,  18,  and  21.  The  curves  on  Figures  12  and  18  were  extracted  from 
page  126  of  Reference  7,  while  the  curves  on  Figures  15  and  21  were 
extracted  from  page  218  of  Reference  7.  For  the  Swerling  I  target 
results  the  curve  could  be  generated  using  Equation  (2.20). 

5 . 3  Probability  of  Detection  Comparison 

The  performance  curves.  Figures  24  through  29,  provide  a  compari¬ 


son  of  the  cell  averaging  and  "greatest-of "  CFAR  probabilities  of 
detection.  They  are  replots  of  Figures  17  through  28  where  the  cell 


Performance  Comparison  of  Cell  Averaging  and  "Greatest 
PFAn  =  TO'4,  Steady  Target,  Square  Law  Detector) 


Figure  26.  Performance  Comparison  of  Cell  Averaging  and  "Greatest-Of " 
(PFAn  =  TO"5,  Steady  Target,  Square  Law  Detector) 


averaging  and  "greatest-of "  CFAR  probabilities  of  detection  curves  are 
combined  for  a  particular  design  probability  of  false  alarm  and  target 
model.  The  actual  probabilities  of  false  alarm  are  given  in  Table  1. 

It  can  be  seen  that  the  cell  averaging  CFAR  has  better  detection 
performance  in  homogeneous  noise  than  the  "greatest-of 11  CFAR.  For  the 
target  model  of  greater  interest,  i.e.,  the  Swerling  I  target,  the 
detection  performance  of  the  two  processors  is  almost  equivalent. 

Since  the  Monte  Carlo  simulation  determines  probabilities  of  detec 
tion,  a  more  meaningful  comparison  could  be  made  using  signal -to-noise 
ratios.  Using  Equation  (3.5),  the  input  signal -to-noi se  ratio  for  the 
cell  averaging  CFAR  and  a  Swerling  I  target  is 


i  -  ffaca-i/n 
-  ' 


(5.6) 


where  N  is  the  number  of  reference  cells,  PFAca  is  the  average  proba¬ 
bility  of  false  alarm,  and  PDcft  is  the  average  probability  of  detection 
Hence,  given  an  average  probability  of  detection,  average  probability 
of  false  alarm,  and  CFAR  window  width,  the  input  signal -to-noise  ratio 
for  the  cell  averaging  CFAR  can  be  determined. 

To  compare  the  two  CFAR  techniques,  an  average  signal -to-noise 
difference  AdB  was  calculated  as  follows.  For  a  particular  PFAp,  i.e., 
1 0-3 ,  10~4,  or  10'5,  CFAR  window  width  N,  i.e.,  8,  16,  or  32,  and 
the  Monte  Carlo  determined  PD' s,  i.e.,  PD^  or  PDg,  the  corresponding 
cell  averaging  CFAR  input  signal -to-noise  ratio  for  both  CFAR  tech¬ 
niques  was  determined  using  Equation  (5.6).  That  is,  the  input 
si gnal -to-noi se  ratio  xrfl (j )  for  the  cell  averaging  is  calculated  as 
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xCA(j) 


l  -  ri/pDrA(j)  PFA n"1/N 


V 


CA _ 

POCA(j)  -  1 


(5.7) 


and  the  input  SNR,  x„(j),  for  the  "greatest-of"  is  calculated  as 


xG(j)  = 


VpMT)  pfad'1/n 


VPDr(j)  -  1 


(5.8) 


where  j  is  an  index  for  the  21  Monte  Carlo  obtained  probabilities  of 
detection  corresponding  to  each  of  the  21  input  signal -to-noise  ratios, 
i.e.,  from  0  dB  to  20  dB  in  increments  of  1  dB. 

After  converting  xCA(j)  and  xG(j)  to  decibels,  an  average  signal- 


to-noise  difference  is  calculated  as 


21  xrA(j)  -  xfi(j) 
AdB  =  £  —  b 

j=l 


21 


(5.9) 


The  results,  given  in  Table  2,  indicate  a  range  of  average  signal- 
to-noise  differences  for  the  Swerling  I  target  as  0.115  dB  to  0.215  dB. 
These  results  are  comparable  to  analytical  results  of  [8]  and  [13]. 


Table  2.  Signal -to-Noise  Ratio  Comparison  of  Cell 
Averaging  and  "Greatest-Of"  CFAR  Processors 


PFAd 

N 

AdB 

8 

BR3SSH 

10"'3 

16 

32 

0.115 

8 

0.215 

10'4 

16 

0.190 

32 

0.142 

8 

10"5 

16 

32 

0.150 
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5.4  Detector  Law  Performance  Comparison 

The  recent  introduction  of  digital  technology  to  radar  signal 
processing  has  necessitated  the  use  of  linear  detectors  for  amplitude 
extraction.  This  is  a  result  of  the  bit  growth  associated  with  a 
squaring  function,  i.e.,  for  B  bit  input  the  output  requires  2-B  bits. 

It  is  further  noted  that  an  exact  linear  detector  is  the  square 
root  of  a  square  law  detector.  Hence,  the  actual  detector  used  is  an 
approximation  to  the  exact  linear  detector.  A  number  of  these  algo¬ 
rithms  have  been  designed  and  normally  take  advantage  of  the  divide  by 
two  which  results  from  right  shifts  of  digital  words.  Two  commonly 
found  algorithms  are 

R  =  MAX  (|I|  ,|Q|)  +  \  MIN  (|I|  JO!)  (5.10) 

and 

R  =  MAX  ( 1 1 1  , 1 Q i )  +  MIN  ( ( 1 1  ,|Q| )  (5.11) 

where  R  is  detector  output  and  |Ij  and  ]Q|  are  the  absolute  values  of 
the  in-phase  and  quadrature  inputs.  Since  there  are  a  number  of 
detector  approximation  algorithms,  no  processor  analysis  is  performed 
using  these  algorithms. 

The  performance  curves  given  in  Figures  30  through  37  were 
obtained  for  the  cell  averaging  and  "greatest-of "  CFAR  processors 
using  exact  square  law  and  linear  detection.  Again  the  probabilities 
of  detection  obtained  by  the  Monte  Carlo  simulation  are  plotted  versus 
the  input  signal -to-noise  ratio.  The  results  obtained  for  a  particular 
CFAR  technique  and  both  detectors  are  plotted  together. 

It  can  be  seen  that  the  performance  of  the  square  law  detector  is 
superior  to  that  of  the  linear  detector  for  any  combination  of  target 


Figure  30.  Performance  Comparison  of  Detector  Laws  for  Cell  Averaging 

(PFAn  =  10~^,  Steady  Target) 


re  31.  Performance  Comparison  of  Detector  Laws  for  Cell  Averaging 
(PFAn  =  10~4,  Steady  Target) 


Figure  32.  Performance  Comparison  of  Detector  Laws  for  "Greatest-Of 

(PFAn  =  TO’3,  Steady  Target) 


Figure  36.  Performance  Comparison  of  Detector  Laws  for  "Greatest 
(PFAn  =  10  Swerling  I  Target) 
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Figure  37.  Performance  Comparison  of  Detector  Laws  for  "Greatest-Of 
(PFAn  =  TO'4,  Swerling  I  Target) 


69 


model,  average  probability  of  false  alarm  and  window  width,  for  either 
CFAR  technique.  For  the  target  of  greatest  interest,  i.e.,  the  Swerling 
I  target,  the  detection  difference  between  using  a  square  law  or  a 
linear  detector  is  small. 

The  actual  PFAs  determined  by  the  Monte  Carlo  simulation  for  the 
linear  detector  system  are  given  in  Table  3. 

Table  3.  Linear  Detector  Probabilities  of  False  Alarms 


PFA0 

N 

PFACA 

PFAg 

10"3 

8 

16 

32 

0.114-2 

0.112-2 

0.98-3 

0.98-3 

0.89-3 

0.89-3 

8 

0.125-3 

0.10-3 

10'4 

16 

0.99-4 

0.93-4 

32 

0.103-3 

0.99-4 

Average  signal -to-noise  ratio  differences  AdB  are  calculated  for  a 
particular  CFAR  procedure  and  the  probabilities  of  detection  obtained 
for  the  two  detectors.  A  positive  AdB  indicates  a  superior  performance 
for  the  square  law  detector.  The  average  signal -to-noise  ratio  differ¬ 
ences  AdB  are  given  in  Table  4  for  a  cell  averaging  CFAR  and  a  Swerling 
I  target  and  in  Table  5  for  a  "greatest-of "  CFAR  and  Swerling  I  target. 
The  results  indicate  that  the  difference  in  the  performance  obtained 
for  either  detector  and  32  reference  cells  is  at  most  0.22  dB.  Hence, 
most  results  obtained  for  a  square  law  system  could  be  used  for  a 
linear  detector  system  as  well.  This  is  an  important  conclusion  since 
the  theoretical  analysis  of  the  CFAR  processors  is  obtainable  only  for 
a  square  law  detector,  and  analysis  costs  are  reduced  by  not  having  to 
simulate  both  detectors. 


'JSL- 


Table  4.  Signal -to-Noise  Ratio  Comparison 
for  Different  Detector  Laws  and 
a  Cell  Averaging  CFAR 


PFAd 

N 

AdB 

10"3 

8 

16 

m 

32 

mm 

10‘4 

M 

Table  5.  Signal -to-Noise  Ratio  Comparison 
for  Different  Detector  Laws  and 
a  "Greatest-of"  CFAR 


PFA0 

N 

AdB 

id-3 

8 

16 

32 

m 

8 

mam 

io'4 

16 

32 

0.21 

5.5  Clutter  Edge  Performance  Comparison 

One  problem  which  must  be  solved  by  adaptive  detection  techniques 
is  tho  regulation  of  false  alarms  in  nonhomogeneous  interference.  For 
radar  mis  would  be  chaff  or  weather  clutter  distributed  in  range. 

The  boundary  of  this  interference,  i.e.,  the  clutter  edge  (Figure  11) 
will  move  into  (or  out  of)  the  reference  cells  as  the  range  cell  of 
interest  approaches  (or  leaves)  the  clutter  area.  Generally,  the 
signal -to-interference  ratio  in  the  clutter  area  is  low  and  the  proba¬ 
bility  of  detection  is  small.  Hence,  the  deviation  of  the  false  alarm 
■  ‘ mr  the  originally  designed  rate  is  of  greater  importance. 
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It  will  be  assumed  that  the  range  extent  of  the  clutter  will  be 
sufficient  to  eventually  cover  all  of  the  CFAR  range  cells,  i.e.,  the 
reference  cells  and  the  cell  of  interest.  The  clutter  will  be  described 
mathematically  as  white  Gaussian  noise  with  the  ratio  of  the  reference 
cell  noise  variance  to  the  cell -of-interest  noise  variance  as 

c2 

Tn  =  n  =  1,  2,  ...  N  .  (5.12) 

a 

It  has  been  shown  [19]  that  for  the  cell  averaging  CFAR 

pfaca  -  n  (l  *  t 


-V 

n  N  ) 


(5.13) 


For  the  condition  where  the  clutter  occupies  <  N/2  reference 
2  2 

cells,  then  rn  =  oJa  =  t  .  In  each  of  these  cells  and  Tn  =  1  for 
non-clutter  cells,  it  follows  that 

™CA  '  ^  ['  +  K'VN  -  1)]  Nl  (5-,4> 


where  PFAp  is  the  homogeneous  interference  design  value. 

For  N/2i.N^<N,  then  the  cell  of  interest  will  also  contain  the 

clutter  and  t  =  1  for  the  clutter  covered  cells.  The  uncovered  cells 
n 

2  2 

will  have  t n  =  a  /a  =  1/t  .  Thus 


(5.15) 


Theoretical  results  for  the  "greatest-of "  CFAR  in  nonhomogeneous  inter¬ 
ference  have  been  performed  only  for  restrictive  cases  [11].  The  Monte 
Carlo  simulation  allows  not  only  for  verification  of  cell  averaging 


theoretical  results  but  also  determination  of  the  "greatest-of" 
performance. 

Performance  comparison  curves  are  shown  in  Figures  38  through  43. 
These  curves  give  the  probabilities  of  false  alarm  versus  the  number  of 
cells  covered  by  the  clutter  for  the  cell  averaging  and  "greatest-of" 
CFAR  processor,,.  The  plots  contain  the  cell  averaging  theoretical 
analysis  and  the  Monte  Carlo  results.  There  are  two  regions  divided  by 
the  cell  of  interest. 

In  region  one,  where  the  clutter  edge  enters  either  CFAR  processor 
window,  the  probability  of  false  alarm  is  smaller  than  the  originally 
designed  probability  of  false  alarm.  Hence,  this  region  is  not 
important  since  both  processors  will  maintain  the  false  alarm  rate 
below  the  design  false  alarm  rate.  For  this  region  the  PFAs  determined 
by  the  simulation  produce  erroneous  or  no  results  for  probabilities  of 

_5 

false  alarm  less  than  10  .  However,  the  cell  averaging  theoretical 

results.  Equation  (5.14),  are  plotted. 

In  region  two,  where  the  clutter  is  in  the  cell  of  interest  and 
there  are  at  least  N/2  +  1  reference  cells,  the  probability  of  false 
alarm  is  now  greater  than  the  design  probability  of  false  alarm  for 
both  processors.  The  worst  case  for  both  processors  is  when  the  cell 
of  interest  and  N/2  +  1  reference  cells  are  covered.  The  actual 
probabilities  of  false  alarm  are  given  in  Table  6. 

It  is  readily  observed  that  both  CFAR  processors  cannot  maintain 
the  design  probability  of  false  alarm  in  certain  clutter  edge  condi¬ 
tions.  However,  the  "greatest-of"  CFAR  is  less  sensitive  to  this 
environment.  From  Table  6  the  probabil ities  of  false  alarm  for  the 
cell  averaging  CFAR  are  a  factor  of  1.4  to  7.4  higher  than  the 
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CA  THEORETICAL 
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NUMBER  OF  CELLS  COVERED  BY  CLUTTER 

Figure  39.  Clutter  Edge  Effects  on  Probability 
of  False  Alarm  (N=8,  tc  =  100,  PFAp  =  10'3) 
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Figure  40.  Clutter  Edge  Effects  on  Probability 
of  False  Alarm  (N=16,  Tc  =  10,  PFAQ  =  10"3) 
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Figure  41.  Clutter  Edge  Effects  on  Probability 
of  False  Alarm  (N=16,  tc  =  100,  PFAq  =  10"3) 
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Figure  42.  Clutter  Edge  Effects  on  Probability 
of  False  Alarm  (N=32 ,  t  =  10,  PFAn  =  10~3) 
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Figure  43.  Clutter  Edge  Effects  Probability 
of  False  Alarm  (N=32,  tc  =  100,  PFAp  =  10'3) 
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"greatest-of"  probabilities  of  false  alarm.  This  factor,  which  shows 
the  difference  in  performance,  is  a  function  of  the  design  probability 
of  false  alarm.  As  shown  in  Moore  and  Lawrence  [11],  for  a  PFAQ  =  10”6, 
N  =  32  and  tc  =  100,  the  cell  averaging  CFAR  probability  of  false  alarm 
is  57.5  times  greater  than  the  "greatest-of"  probability  of  false  alarm. 
This  feature  is  the  primary  advantage  that  the  "greatest-of"  CFAR  has 
over  the  cell  averaging  CFAR. 

Table  6.  Comparison  of  Clutter  Edge  Probabilities 
of  False  Alarm  (PFAq  =  10"3) 


N 

pfaca 

pfag 

8 

■ 

.93-2 

0.66-2 

16 

0 

.15-1 

0.54-2 

32 

0 

.18-1 

0.35-2 
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uantization  Consideration 


Recent  radar  signal  processors  are  implemented  digitally,  hence, 
the  effect  of  quantization  noise  must  be  considered  when  specifying  a 
desired  probability  of  false  alarm.  For  this  analysis  the  quantization 
or  A/D  conversion  will  occur  after  the  detector  with  wordlengths  of  6, 
8,  and  10  bits  and  the  wordlength  is  not  truncated  in  the  CFAR 
processor.  Since  a  given  CFAR  processor  could  have  any  combination  of 
wordlengths,  truncation  schemes  and  assumed  saturation  level,  no  effort 
will  be  made  to  determine  a  general  method  to  maintain  a  given  false 
alarm. 

The  quantization  errors  for  8  and  10  bits  wordlengths  are  negli¬ 
gible,  thus,  only  the  results  for  6  bits  are  given  in  Tables  7 


FkZ 


through  10.  Tables  7  and  8  give  the  probabilities  of  false  alarm  for 

2  2 

the  two  CFAR  processors  in  homogeneous  noise  with  a  =1  and  a  =2  and 
6  bits  of  quantization.  Table  9  gives  the  average  signal -to-noise 
difference,  based  on  the  cell  averaging  CFAR  as  described  in  this 
chapter,  for  both  CFAR  processors  and  6-bit  quantization.  Table  10 
gives  probabilities  of  false  alarm  for  both  processors,  linear  detector, 
and  6  bits. 


Table  7.  Quantization  Effects  on  Probability 

o 

of  False  Alarm  (Square  Law,  6  bits,  a  =1) 


PFAD 

N 

PFACA 

pfag 

8 

0.55-3 

0.52-3 

10  3 

16 

0.57-3 

0.71-3 

32 

0.78-3 

0.82-3 

8 

0.19-4 

0.51-4 

10  4 

16 

0.51-4 

0.57-4 

32 

0.64-4 

0.66-4 

c 

8 

_ 

0.27-5 

10~5 

16 

0.19-5 

0.42-5 

32 

0.36-5 

0.64-5 

Table  8.  Quantization  Effects  on  Probability 

2 

of  False  Alarm  (Square  Law;  6  bits,  o  =  2) 


pfad 

N 

pfaca 

pfag 

8 

0.95-3 

0.850-3 

10'3 

16 

0.94-3 

0.103-2 

32 

0.10-2 

0.930-3 

A 

8 

0.68-4 

0.910-4 

10'4 

16 

0.93-4 

0.105-3 

32 

0.106-5 

0.112-3 

c 

8 

0.50-5 

0.670-5 

10'5 

16 

0.77-5 

0.860-5 

32 

0.79-5 

0.108-4 

Table  9.  Quantization  Effects  on  Signal -to-Noise 
Ratio  (Square  Law,  6  bits) 
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PFAd 

N 

aca 

ag 

o 

8 

0.362 

0.418 

10'3 

16 

0.286 

0.116 

32 

0.264 

0.248 

8 

0.460 

0.256 

10‘4 

16 

0.314 

0.257 

32 

C.280 

0.217 

8 

0.600 

0.382 

10"5 

16 

0.365 

0.276 

32 

0.294 

0.240 

Table  10.  Quantization  Effects  on  Probability  of 

2 

False  Alarm  (Linear  Detector,  6  bits  o  =1) 


pfad 

N 

pfaca 

PFAg 

8 

0.117-2 

0.75-3 

io'-5 

16 

0.107-2 

0.89-3 

32 

0.980-3 

0.86-3 

8 

0.124-3 

0.101-3 

10"4 

16 

0.108-3 

0.970-4 

32 

0.108-3 

0.100-3 

It  is  observed  from  Table  7  that  the  probability  of  false  alarm 
has  decreased  for  both  processors  which  is  undesirable  since  this  is 
not  the  original  design  value.  Even  if  the  CFAR  scale  factors  K  and  Kg 
are  adjusted  to  obtain  the  design  probability  of  false  alarm,  the  proba¬ 
bility  of  false  alarm  would  change  if  the  standard  deviation  of  the 
noise  changed  as  shown  in  Tab'le  8.  Since  the  probability  of  false 
alarm  (Table  7)  decreased  due  to  quantization,  the  probability  of 
detection  is  also  reduced.  This  is  s  own  in  Table  9  where  the  average 
signal -to-noi se  differences  between  no  quantization  and  quantization 
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are  given  for  both  processors.  Table  10  shows  no  real  difference  in 
probability  of  false  alarm  due  to  quantization  because  of  the  reduced 
dynamic  range  of  the  linear  detector  outputs. 

The  effects  of  A/D  quantization  and  finite  wordlengths  must  be 
considered  when  implementing  a  digital  CFAR  processor.  Even  though  some 
analytical  effort  has  been  performed  [20],  a  more  complete  study  can  be 
achieved  only  through  simulation. 

5.7  Non-Gaussian  Interference  Results 

The  cell  averaging  CFAR  and  "greatest-of"  CFAR  processors  assume 
the  noise  amplitude  distribution  is  Gaussian  with  an  unknown  power. 

In  several  instances  this  is  not  a  valid  assumption  and  a  changing 
probability  density  function  can  be  encountered  due  to  a  lack  of 
clutter  rejection  by  the  MTI. 

Several  investigations  of  natural  clutter  characteristics  have 
shown  that  clutter  returns  can  be  described  by  log-normal  or  Wei  bull 
[18]  types  of  distributions  where  the  Weibull  pdf  includes  the 
Rayleigh  pdf  as  a  special  case. 

The  Weibull  pdf  is  a  single  variate  function  having  two 
parameters,  a  and  b,  and  is  given  by 


where  o°  is  the  variate  in  terms  of  the  clutter  backscatter  coefficient, 
b  =  1/A  (A  =  Weibull  slope  parameter)  and 


where  a0  =  median  value  of  Weibull  pdf. 
m 


(5.17) 
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The  probability  of  false  alarm  performance  against  Weibull  for 
A  =  1,  2,  3  for  both  CFAR  processors  is  given  in  Tables  11  through  13. 

For  A=1 ,  the  Weibull  pdf  reduces  to  the  exponential  pdf,  therefore, 
the  probability  of  false  alarms  obtained  are  the  originally  designed 
values . 

For  A=2  and  A=3,  which  are  representative  of  natural  clutter  [18], 
the  probabilities  of  false  alarm  increase  by  a  factor  of  approximately 
100  and  1000,  respectively,  for  both  processors.  This  increase  is 
unacceptable. 


Table  11.  Probabilities  of  False  Alarm  in  Presence 
of  Weibull  Clutter  (A=l ) 


PFAd 

N 

pfaca 

pfag 

10-3 

8 

0.112-2 

0.890-3 

16 

0.117-2 

0.118-2 

32 

0.112-2 

0.106-2 

8 

0.119-3 

0.125-3 

io'4 

16 

0.106-3 

0.111-3 

32 

0.110-3 

0.129-3 

8 

0.102-4 

0.75-5 

10'5 

16 

0.133-4 

0.142-4 

32 

0.117-4 

0.118-4 

Table  12.  Probabilities  of  False  Alarm  in  Presence 
of  Weibull  Clutter  (A=2) 


pfad 

N 

pfaca 

pfag 

10'3 

8 

16 

0.331-1 

0.298-1 

0.252-1 

0.250-1 

32 

0.283-1 

0.237-1 

8 

0.170-1 

0.136-1 

10‘4 

16 

0.164-1 

0.129-1 

32 

0.154-1 

0.125-1 

io'5 

8 

16 

0.878-2 

0.916-2 

0.686-2 

0.699-2 

32 

0.905-2 

0.717-2 
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Table  13.  Probabilities  of  False  Alarm  in  Presence 
of  Weibull  Clutter  (A=3) 
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PFAd 

N 

mQ 

8 

0.619-1 

0.500-1 

10"3 

16 

0.498-1 

0.410-1 

32 

0.435-1 

0.356-1 

A 

8 

0.431-1 

0.355-1 

10'4 

16 

0.359-1 

0.291-1 

32 

0.307-1 

0.248-1 

c 

8 

0.304-1 

0.246-1 

10"5 

16 

0.265-1 

0.210-1 

32 

0.229-1 

0.181-1 

A  CFAR  processor  has  been  designed  which  maintains  false  alarm 
regulation  in  log-normal  and  Weibull  clutter  [21].  Also,  a  Weibull 
loss  has  been  presented  for  the  cell  averaging  CFAR  designed  to  main¬ 
tain  a  constant  false  alarm  rate  in  various  Weibull  clutter  environ¬ 
ments  [16]. 

5.8  Interfering  Target  Results 

The  detection  performance  of  both  CFAR  processors  will  be  affected 
by  a  target  or  targets  occupying  the  CFAR  window  when  a  target  is  in 
the  cell  of  interest.  The  interfering  target(s)  can  reduce  the  proba¬ 
bility  of  detection  to  an  unacceptable  value  as  shown  by  Finn  and 
Johnson  [14]  in  their  Figure  18  for  a  square  law  detector  and  a  target 
pair  using  a  cell  averaging  CFAR. 

For  a  limited  detection  performance  comparison,  two  Swerling  I 
targets  will  be  assumed:  white  Gaussian  noise  and  a  square  law 
detector.  Three  cases  will  be  simulated:  1)  the  target  of  interest 
is  10  dB  above  the  noise  and  the  interfering  target  is  7  dB  above  the 
noise,  2)  the  target  of  interest  and  the  interfering  target  are  both 
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10  dB  above  the  noise,  and  3)  the  target  of  interest  is  10  dB  above  the 
noise  and  the  interfering  target  is  13  dB  above  the  noise. 

The  probabilities  of  detection,  and  PEFg,  determined  by  the 
simulation  are  given  in  Tables  14  through  16.  Using  Equation  (3.5)  and 
the  probabilities  of  detection,  IF  signal -to-noise  ratios,  SNRCA  for 
the  cell  averaging  CFAR  and  SNR^  for  the  "greatest-of"  CFAR  are 
calculated.  The  last  column,  A^NR,  provides  a  measure  of  performance 
comparison  between  the  two  processors. 

It  is  readily  observed  that  the  detection  performance  decreases  as 
N  decreases.  Even  for  a  large  N,  i.e.,  N=32,  the  cell  averaging  CFAR 
suffers  a  detectability  loss  of  0.7  dB,  1.2  dB,  and  2.0  dB  for  cases 
1,  2,  and  3,  respectively.  The  signal -to-noise  difference,  ASNR,  gives 
the  amount  the  input  SNR  could  be  reduced  for  the  cell  averaging  CFAR 
and  maintain  equivalent  performance  with  the  "greatest-of"  CFAR.  The 
range  of  A<-NR  is  from  0.3  dB  to  0.7  dB. 

While  both  processors  are  sensitive  to  an  interfering  target 
environment,  the  cell  averaging  CFAR  is  superior  to  the  "greatest-of" 
in  this  type  of  environment.  This  advantage  would  have  to  be  consid¬ 
ered  when  designing  a  CFAR  processor. 


Table  14.  Signal -to-Noise  Ratio  Comparison 
for  a  7  dB  Interfering  Target 


PFA0 

N 

^CA 

snrca 

SNRg 

asnr 

8 

0.266 

8.2 

0.239 

7.8 

0.4 

10'3 

16 

0.381 

8.9 

0.357 

8.5 

0.4 

32 

0.449 

9.3 

0.431 

9.0 

0.3 

8 

0.140 

8.3 

0.124 

8.0 

0.3 

10~4 

16 

0.259 

8.9 

0.236 

8.6 

0.3 

32 

0.334 

9.4 

0.317 

9.1 

0.3 

8 

0.066 

8.4 

0.058 

8.1 

0.3 

10'5 

16 

0.166 

9.0 

0.148 

8.6 

0.4 

32 

0.246 

9.4 

0.230 

9.1 

0.3 

Table  15.  $ignal-to-Noise  Ratio  Comparison 
for  a  10  dB  Interfering  Target 
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PFA0 

N 

P0CA 

SNP 

CA 

.  <-3 

IS 

snrg 

asnr 

8 

0.194 

7.0 

0.155 

6.3 

0.7 

10'J 

16 

0.318 

8.0 

0.278 

7.4 

0.6 

32 

0.406 

8.7 

0.374 

8.3 

0.4 

8 

0.097 

7.3 

0.081 

6.8 

0.4 

10'4 

16 

0.205 

8.1 

0.175 

7.6 

0.5 

32 

0.294 

8.8 

0.264 

8.4 

0.4 

c 

8 

0.044 

6.4 

0.036 

6.1 

0.3 

io"5 

16 

0.126 

8.2 

0.106 

7.8 

0.4 

32 

0.210 

8.8 

0.183 

8.4 

0.4 

Table  16.  Signal -to-Noise  Ratio  Comparison 
for  a  13  dB  Interfering  Target 


™D 

N 

PDCA 

SNRCA 

snrg 

asnr 

8 

0.124 

5.6 

0.099 

5.0 

0.6 

10'3 

16 

0.243 

6.9 

0.196 

6.1 

0.8 

32 

0.348 

7.9 

0.299 

7.2 

0.7 

8 

0.056 

3.5 

0.046 

2.8 

0.7 

10"4 

16 

0.146 

5.1 

0.116 

4.4 

0.7 

32 

0.242 

6.3 

0.201 

5.7 

0.6 

8 

0.024 

6.4 

0.019 

6.1 

0.3 

10~5 

16 

0.085 

7.2 

0.066 

6.2 

0.5 

32 

0.164 

8.1 

0.133 

7.5 

0.6 

5.9  Summary 

The  detection  performances  of  the  cell  averaging  and  "greatest-of" 
CFAR  have  been  determined  and  presented.  Both  CFAR  processors  were 
designed  for  three  common  probabilities  of  false  alarm  and  the  actual 
probabilities  of  false  alarm  were  obtained  by  the  Monte  Carlo  simula¬ 
tion.  For  a  given  probability  of  false  alarm,  the  two  processors  have 
essentially  equivalent  detection  performance  with  the  "greatest-of " 
having  approximately  0.2  dB  loss  as  compared  to  the  cell  averaging. 
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While  both  processors  have  only  a  slight  degradation  when  a  linear 
detector  instead  of  a  square  law  detector  is  used,  their  performance  is 
unacceptable  in  Weibull  clutter  and  is  affected  by  finite  wordlength 
processing. 

The  two  main  areas  of  performance  comparison  are  the  probability 
of  false  alarm  regulation  in  clutter  edges  and  the  probability  of  detec¬ 
tion  in  an  interfering  target  situation.  The  "greatest-of"  proved  to 
regulate  false  alarms  much  better  in  the  clutter  edge  environment  while 
causing  an  additional  detection  loss,  between  0.3  and  0.7  dB,  in  the 
interfering  target  environment  as  compared  to  the  cell  averaging  method. 

It  is  obvious  that  the  selection  of  either  the  cell  averaging  CFAR 
or  the  "greatest-of"  CFAR  should  be  based  on  the  expected  radar  envir¬ 
onment.  Due  to  similarity  in  their  implementation,  a  combination  of 
the  two  processors  and  supporting  selection  logic  could  provide  an 
overall  improved  CFAR  performance. 


CHAPTER  VI.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 
6.0  Summary 

A  computer  simulation  has  been  developed  for  a  performance  compari 
son  of  two  commonly  known  CFAR  techniques.  The  two  techniques  are  the 
cell  averaging  and  "greatest-of. "  The  comparison  is  based  on  the  proba 
bili ties  of  detection  and  the  probabilities  of  false  alarm  obtained  by 
performing  Monte  Carlo  passes  of  the  simulation. 

The  two  CFAR  processors  were  designed  for  average  probabilities 
-3  -4  -5 

of  false  alarm  of  10  ,10  ,  and  10  .  These  false  alarm  rates  were 

verified  by  the  simulation.  Probability  of  detection  versus  input 
signal-to-noise  ratio  curves  for  each  false  alarm  rate  were  generated 
for  both  processors.  Two  target  models  were  used:  the  steady  or  non¬ 
fluctuating  target  and  the  Swerling  I  target.  The  probability  of 
detection  results  were  utilized  to  make  a  signal-to-noise  ratio 
difference  comparison  and  indicated  that  the  cell  averaging  CFAR  would 
require  approximately  0.2  dB  less  input  signal-to-noise  ratio  for 
equivalent  performance  to  the  "greatest-of"  CFAR. 

The  detection  performance  for  a  linear  detector  system  was  deter¬ 
mined  for  both  CFAR  processors  and  target  models.  Comparing  the 
results  to  a  square  law  detector  system  indicates  that  the  detection 
performance  is  degraded  by  the  use  of  a  linear  detector.  However, 
this  degradation  is  negligible,  especially  for  a  Swerling  I  target. 

Two  important  analyses  were  performed:  thn  clutter  edge  per¬ 
formance  comparison  and  the  interfering  target  results. 
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A  clutter  edge,  i.e.,  residual  clutter  distributed  in  range, 
affects  both  CFAR  processors  by  increasing  the  probability  of  false 
alarm  above  the  design  value.  The  amount  the  false  alarm  rate 
increases  is  a  function  of  the  design  false  alarm  rate  and  the  CFAR 
window  size.  The  "greatest-of"  technique  provides  better  false  alarm 
control  than  the  cell  averaging  technique  in  a  clutter  edge  condition. 

An  interfering  target,  i.e.,  a  target  which  is  in  a  reference  cell 
of  the  CFAR  processor,  degrades  the  detection  performance  of  both  CFAR 
processors.  In  general,  the  amount  of  degradation  is  a  function  of 
the  interfering  target  power,  the  design  probability  of  false  alarm  and 
the  CFAR  window  size.  The  cell  averaging  technique  provides  better 
detection  performance  than  the  "greatest-of"  technique  in  an  interfer¬ 
ing  target  environment. 

The  quantization  analysis  demonstrated  that  the  probability  of 
false  alarm  and  the  probability  of  detection  are  affected  by  finite 
wordlength  arithmetic.  For  some  finite  wordlength  CFAR  processors,  the 
false  a’ arm  rate  will  change  as  the  interference  power  changes. 

The  non-Gaussian  interference  analysis  demonstrated  the  unaccept¬ 
able  false  alarm  rates  obtained  in  Weibull  interference  for  both  CFAR 
processors. 

6.1  Conclusions 

The  cell  averaging  and  "greatest-of"  CFAR  processors  can  be 
designed  to  maintain  a  constant  false  alarm  rate  in  homogeneous  white 
Gaussian  noise.  The  probability  of  detection  obtained  for  a  given 
probability  of  false  alarm  increases  as  the  number  of  CFAR  processor 
reference  cells  increase  for  both  processors.  The  two  CFAR 
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techniques  have  almost  equivalent  performance  in  white  Gaussian  noise 
with  the  cell  averaging  CFAR  having  a  slight  advantage. 
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The  simulation  results  have  shown  a  negligible  improvement  obtained 
for  a  square  law  detector  over  a  linear  detector.  Hence,  the  analytical 
results  developed  for  the  two  CFAR  processors  and  a  square  law  detector 
could  be  used  to  describe  the  performance  for  a  linear  detector  system. 

Whereas  both  processors  fail  to  maintain  the  design  probability 
of  false  alarm  in  a  clutter  edge  environment,  the  "greatest-of "  tech¬ 
nique  is  affected  less  than  the  cell  averaging  technique  and  should  be 
a  prime  CFAR  candidate  if  such  an  environment  is  anticipated. 

An  interfering  target  will  degrade  the  detection  performance  of 
both  CFAR  processors.  The  cell  averaging  technique  is  affected  less 
than  the  "greatest-of"  technique  and  should  be  a  prime  CFAR  candidate 
if  interfering  targets  are  considered  to  be  a  dominant  problem. 

Finally,  the  performance  of  both  processors  is  affected  by  finite 
wordlength  arithmetic  and  the  phenomenon  should  be  analyzed  when  imple¬ 
menting  either  CFAR  technique.  The  unacceptable  false  alarm  rates 
obtained  for  both  CFAR  processors  when  Weibull  clutter  is  in  the 
reference  cells  requires  utilization  of  a  different  CFAR  if  this  is  the 
expected  environment. 

The  results  agree  with  those  previously  available.  But  the  inter¬ 
fering  target  performance  comparison  and  the  "greatest-of"  performance 
in  Weibull  clutter  and  non-Gaussian  interference  represent  results 
presently  not  available. 
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6.2  Recommendations 

The  recommendations  for  future  work  are  to: 

1)  Improve  the  random  number  generator  so  that  probabilities 

_5 

of  false  alarms  less  than  10  can  be  verified. 

2)  Use  measured  radar  data  as  an  input  to  the  simulation  to 
compare  the  CFAR  processors. 

3)  Perform  an  extensive  study  of  linear  detector  approxima¬ 
tion  algorithm's  effect  on  CFAR  performance. 

4)  Perform  an  extensive  study  of  finite  wordlength  effect 
on  design  of  and  performance  of  CFAR  processors. 

5)  Determine  realistic  models  of  jammers  and  perform  a 

study. 

6)  Develop  environmental  models  which  contain  clutter  edges 
and  interfering  targets  and  determine  selection  logic  for  the 
"greatest-of"  and  cell  averaging  CFAR  processors  to  optimize  CFAR 
performance.  Utilization  of  tracking  information  should  be  considered. 
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APPENDIX  A 

FIXED  THRESHOLD  PROBABILITY  DENSITY  FUNCTIONS  DERIVATIONS 


A.  1  Introduction 

This  appendix  gives  derivations  for  the  probability  density  func¬ 
tions  for  the  square  law  or  linear  detected  steady  target  and  for  the 
square  law  detected  Swerling  I  target.  These  probability  density 
functions  are  common  equations  and  the  derivations  can  also  be  found 
in  Marcum  [1]  for  the  steady  target  and  Swerling  [2]  for  the  Swerling 
I  target. 

The  characteristic  function  approach  is  used  to  obtain  the  pdf  for 
the  square  law  detected  output.  The  steady  target  is  assumed  to  be 
distributed  in  the  I  and  Q  channels  by 

Sj  =  P  cos  (e) 

Sq  =  P  sin  (e)  (A.  1 ) 

and 

y  -  (Sj  +  Xj)  +  (Sg  +  Xq) 

and  this  changes  the  pdf's  for  Xj  and  Xq  by  a  shift  to  these  mean 
values,  i.e.,  zero  mean  Gaussian 


p(xj  =  — —  exp 
y^rra 


2a2 


2-, 


I,  Q 


(A. 2) 


In  the  following  derivation  the  subscript  k  will  be  dropped  on  x^ 
and  S^  to  simplify  the  equations. 

2 

The  characteristic  function  for  x  is 


! 
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<j>  ?(v)  =  E{exp(jvx2)}  =  — L-  f 
%c  V2tto  j~co 


exp  [-  (»-S)Z-  Jv»Vl  dx 

2a 


f"  exp [-  dx  .  (A. 3) 

I  2  a 


2 

Now  let  C  =  (1  -  jv2a  )  and  complete  the  square. 


*  2<*>  '  -j=~  I 

x  y2ir  a 


2  2S  .  si  si  +  s£' 

x  “  T  C  “  r2  r2 
exp - 5 - - - —  dx 


(2a  /C) 


S2  S2 


-4  exp  5^-1  — J-5-  f  expf-  -t^l  dx  . 
VC  |_  (2a/C)  J  y^TT  L  (2a  /C) 

VC  (A. 4) 


After  integration. 


<P  P(v)  =  ~  exp 


Consequently,  the  characteristic  function  for  y  is 


4>  (v)  =  <J>  2(v)  4>  2(v)  =  £  exp 


- 1  (si  +  so)(r-  0 


XI  XQ 


iexp[-^]exp[^_  = 


{’SK 


P2  1 
2?  (1  -  j2a2vT- 


1  - j2a  v 


(A. 5) 


Let  x  =  P2/(2a2)  and  -v  =  u,  then 


exp(-x) 
2tt  ( 2  a2 ) 


T  x  1 

exp  —j- - T 

[2 alju  +  l/2a  . 


ju  +  — y 

2a 


exp[juy]du 


-  expf -4,1  I 


2  o  ' 
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Io(2”f)u(y>  (A-6> 

This  pdf  is  the  well  known  Rician  distribution  [3].  A  linear  detector 
is  given  as  z  =  /y  and  by  a  change  in  variables  the  pdf  for  a  linear 
detector  can  be  obtained  from  Equation  (A. 6). 


2o 


exp 


A. 2  Derivation  of  Probability  Density  Function  for  a 
Square  Law  Detected  Swerling  I  Target  Plus  Noise 

When  a  Swerling  I  target  is  assumed,  the  results  for  the  pdf  are 

simplified.  This  tvpe  of  target  assumed  that  the  group  of  N  returns 

have  a  constant  signal -to-noise  ratio  but  that  from  group  to  group  the 

pdf  is 

p( x )  =  ~  exp[-x/ x  ]u(x)  (A. 7) 

x 


O  O  _ 

where  x  =  P  /(2a  )  and  x  is  the  average  signal -to-noise  ratio.  This  is 
used  with  the  characteristic  function  of  Equation  (A. 5)  to  obtain 


\(v) 

Vv) 


Vv) 


oo 

=  E{4>y(v)}  =  j  4>y(v)p(x)dx 

=  /0  T  exPf"xJ  exP  [f]  exP 

=  —  f  exp  x  (-1  +  J-  -  ~) 
cx  J0  L  v  u  X  / 


dx 


dx 


cx(^7'') 


[exp(-°°)  -  exp(O)] 


-1 


X  -  c  -  cx 


-1 


x  -  1  +  j2a2  v  -  x  +  j2a^vx  1  -  jZ a2  v(l  +  x) 


1 


1 


- - - -  T 

2  a  ( 1  +  x )  -jv  +  — p 

n  > 


2a  (1  +  x) 


(A. 8) 


The  pdf  for  y  can  be  determined  from  this  to  be 
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APPENDIX  B 

DERIVATION  OF  THE  GREATEST-OF  CFAR  PERFORMANCE  EQUATIONS 


B.O  Introduction 

This  appendix  gives  the  derivations  for  the  probability  of  false 
alarm  and  probability  of  detection  equations  for  the  "greatest  of" 

CFAR.  The  derivations  were  originally  derived  by  Moore  [8],  It  should 
be  noted  that  in  an  independent  concurrent  effort,  Hansen  and  Sawyer 
[13]  have  derived  the  same  equations. 

B . 1  Derivations  of  the  Greatest-Of  CFAR  Performance  Equations 

In  the  "greatest-of"  CFAR  method  two  independent  thresholds  are 
calculated,  then  the  largest  one  is  selected,  viz., 


M 


Y,  = 


’  H  A 


Y2  =  «  j?,  yj 


i=l 

M 


Yth  =  MAX  (Yr  Y2)  (B.l) 

where  a  simplified  notation  is  used  for  the  subscripts  on  y  and  it  is 
implied  that  the  ranges  of  the  summations  are  M  but  that  i=j. 

The  pdf  descriptions  of  Y-j  and  Y^  can  be  given  by 


P(Y)  •(») 


2  o 


1 


2  uny 


r&r 


exp 


MY 

K2?J 


u(Y) 


(B.2) 


Papoulis  [9]  gives  an  expression  for  finding  a  pdf  of  the  maximum  of 
two  random  variables,  cf . ,  Equation  (15),  p.  193, 

PG0<V  '  2Fmp(Y)l  -  ?FY(Yth)PY(Yth) 

th 

where  Fy(-)  is  the  cumulative  distribution  unction  for  Y. 


(B.3) 
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This  could  be  used  with  fixed  threshold  probability  of  false  alarm 
Equation  (2.8)  to  obtain  the  expected  PFA,  i.e.. 


PFA 


GO 


CO 

f  exp 

1 

1 

-< 

J 

J0 

L  2a"  J 

PG0<¥th>d¥th 


=  2  f  exp 

J0 


-V 


th 


2a 


Fv<¥th)Fv'Yth>dYth 


This  can  be  integrated  by  parts  as  follows: 

u  =  Fv(Vth)  to  •  F;(Yth)dVth 


dv  =  exp 
v  -f, 


r-y 


th 


L  2o 

Y 


exp 


2 

th 


2a 


F;<¥th>dvth 

»y<¥th>d¥th 


(8.4) 


(B.5) 


Note  that  v  is  in  the  same  form  used  to  obtain  PFA  for  a  cell  averaging 
CFAR  Equation  (3.3),  but  this  is  not  a  definite  integral.  It  follows 
that 


(B.6) 


where  PFA  is  the  value  associated  with  a  conventional  cell  averaging 
CFAR  of  window  size  M.  Thus 
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pfago  -  *K  '-H  ‘  2FV<Vvl 

♦  v  HT  r/ithV 


'0  \  2a 


0  +  k)  py^Yth^dYi 


The  first  term,  FY(Vth)v,  yields  0,  thus 


PFAgq  =  2PFA 


/  M\m 

L  -TH--  X  2c'2'*  exp[-a(,  +  ^) 


Py(2o  a)da 


/mxM  ,  M-l  (l  +  JT) 

■»®  bAtj yi 


m+M-1 

a  exp 


[-•('♦?)} 


__/mxM  _ 2 _ M-l  (m+M-  1)! 

=  PFA(")  * f\' 

_  2PFA  Mf 1  (m  +  M  -  1)1 

/  i/\M  m_n  , 

(M  -  1 ) !  (2  +  ^)  m-°  /2  +  J\ 

mil  - rr  1 

\1  +  -/ 

V  +  M' 


/  K\  _ '1/M 

Since  (l  +  ^ j  =  PFA  ,  then 


(B.7) 


pp7C  _  2PFA _ (m  +  M-1)! 

GO  ”  ,  /  _ -1/MvM  %  f  _ -1  /M\ m" 

(M  -  1 )!  (l  -  PFA  )  m  U  m!  (l  +  PFA  ) 

This  relates  the  probabilities  of  false  alarm. 

The  probability  of  detection  for  a  Swerling  I  target  with 
"greatest-of"  CFAR  is  given  by 


(B.8) 


<*Vv 
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co  r- 

1,60  =  fQ  exp[' 


’th 


2a"  (1  +x)J 


2  f  2a  exp 

J0 


a 

L  1  +  x 


PG0<VdVth 


Fy(2a^a)  Py(2a^a)da 


Integration  by  parts  gives 
2 


u  =  Fy(2a  a)  du  *  Fy (2a^a)  2a  aa 
dv  «  2a"  exp^-1~Jpy(20^a)da 

/M\M  1  „  f  /M  ,  1  \ 

v  “  (k)  TfTTJT  J0a  exp['U  +  1  +-j 


da 


-  ® 


1 

Tm^TTF 


V*. 

i  r 

(u 

|M  Jo 

\  M(1  +  x)> 

TirrrjT  /  \m  exp(-b)  J0  m!  • 


'1  + 


M(  7 


K_\" 

+  x)/ 


where 


(B.9) 


(B.10) 


(B.l 1 ) 


Thus,  using  the  probability  of  detection  for  an  equivalent  sized  cell 
averaging  CFAR, 


M-l  Km 

V  =  -P„  exp(-b)  I  2t 

m=u 


(B.12) 


and 


9  —  9  /•  h  9 

PD,G0  '  2ry(2o  a)v  ^  -  2Pd  £  ^  fg  «P<-»>  PV<2°  a)da 

H_!  /M  +  _J\ 

2PD  fM\M  I  VK  lii/  a/2M  +  _jJ 

-  VK  l  +  * 


(M^TTI  Vk/  n=0  nil 


C  m+M-l 
I  a  exp 

J0 


da 


D,GO  (M  -  TT1  \K 


+  x/ 


m=0 


nt! 


<?  •  *) 


M+m 


I 


m+H- 1 


exp[-c]dc 


2Pr 


D 

Im^DT 


1  Mf]  (m+M-1)! 

vM  ** 


2  + 


— ) 
+  X  )  ' 


m=0 


M(  1  +  x) 


m! 


2  + 


K 


M  ( 1  +  x) 


1  + 


K 


M(1  +  x )  J 


2P 

0 

TTTTyr 


i 


M-l 


(m+M-1)! 


(1  +  Pd‘1/M)M  m~°  m!  (1  +  Pd'1/M) 


(B. 


This  is  the  desired  expression  for  the  probability  of  detection. 
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APPENDIX  C 

MONTE  CARLO  RUN  ESTIMATION 

The  utilization  of  Monte  Carlo  simulations  for  estimation  of 
probabilities  of  false  alarm  and  probabilities  of  detection  has  a 
statistical  uncertainty  associated  with  it.  This  appendix  determines 
values  required  to  give  a  priori  probabilities  for  a  specified  range 
of  the  estimated  parameter. 

Let  yn  represent  the  nth  target/no  target  decision  for  the  cell 
of  interest.  Thus  yp  will  equal  either  0  or  1 .  The  probability  that 
yn  =  1  will  be  denoted  as  p.  Two  cases  are  considered,  viz.,  noise- 
orily  and  signal -plus-noise. 

Thus 

Prob  [yp  =  1  |  noise  only]  =  p  =  PFA 

Prob  [y  =  1  |  signal -plus-noise]  =  p  =  PD  (C.l) 

an  estimate  of  p  can  be  formed  by  calculating  the  arithmetic  mean  of 

N  determinations,  i.e., 

*  'IT  Vn  '  (C'2> 

n=  l 

This  represents  an  unbiased,  efficient  and  consistent  estimator  to  the 
expected  value  of  y  .  It  is  possible  to  obtain  the  mean  and  variance 
of  y  in  terms  of  the  mean  and  variance  of  y  ,  i.e., 

Ely)  =  J-  I  E{yn)  =  E{yn)  =  p 

E|(y-p)2)  =  VAR(y)  =  -AR^-y-n-)  =  •  (C.3) 
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One  approach  to  finding  the  required  value  for  N  is  to  use  Chebyshev's 
Inequal ity,  i  .e. , 


P[|y  -  E|y)|  >  e]  il 


(C.4) 


letting  the  value  of  e  depend  on  p,  i.e.. 


e  =  kp 


(C.5) 


yields 


p[-e  <  y  -  p  <  e]  >  1  -  --  —  = 


(C.6) 


Thus  if  the  estimate  (y)  to  p  is  to  be  within  some  specified  range  of 
p(±e)  with  better  than  some  specified  probability  (K)  then  Equation 
(C.6)  can  be  used  to  determine  the  sufficient  value  for  N.  Typical 
results  are  given  in  Table  C.l. 

Table  C.l.  Values  of  N  Obtained  by  Chebyshev's  Inequality 


K  =  0.5 


K  =  0.9 


0.6  1 .33  x  10 


0.7  8.57  x  10' 

0.8  5  x  103 

0.9  2.22  x  10' 


1  x  10 
6.67  x  104 
4.29  x  104 
2.25  x  10^ 
1.11  x  104 


1  x  10 
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APPENDIX  D 

FLOATING  POINT  SYSTEMS  AP-120B 


The  Floating  Point  Systems  AP-120B  is  a  loosely  coupled  synchronous 
array  processor  which  uses  pipelined  arithmetic  elements.  The  array 
processor  uses  a  38-bit  floating-point  format  and  has  a  cycle  time  of 
167  nsec.  Figure  D.l  shows  the  structure  of  the  AP-120B,  which  con¬ 
sists  of  an  interface  to  the  host  computer,  a  program  memory,  a  16-bit 
integer  ALU,  data  memory,  table  memory,  accumulators,  I/O  interface, 
and  arithmetic  elements. 
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The  interface  controls  data,  program  transfer,  and  format  conver¬ 
sion  between  the  host  and  array  processor. 

Control  consists  of  two  elements:  a  16-bit  ALU  which  performs 
integer  address  indexing  and  loop  counting  for  all  of  the  memory 
elements.  The  second  element  is  the  program  memory  which  contains  the 
microcode  to  be  executed  in  the  array  processor.  This  memory  is  64 
bits  wide  with  each  word  being  subdivided  into  10  command  fields.  Each 
command  field  controls  an  element  in  the  array  processor,  thus  every 
element  can  be  active  in  every  machine  cycle. 

The  main  data  memory  is  used  for  data;  the  table  memory  is  used 
for  storing  constants  and  the  accumulator  blocks  for  intermediate 
result  storage. 

The  AP-120B  uses  a  unique  bus  structure  in  that  there  are  dedi¬ 
cated  paths  between  each  memory  and  each  arithmetic  element,  thus  maxi¬ 
mizing  the  flow  of  operands  and  resultants  between  functional  elements. 

The  I/O  interface  allows  the  attachment  of  peripheral  devices 
directly  to  the  array  processor. 

The  arithmetic  elements  consist  of  a  3-stage  multiplier  and  a 
2-stage  adder,  each  stage  running  at  the  cycle  time  of  the  array 
processor  (167  nsec),  thus  a  multiply-add  can  be  obtained  in  every 
cycle  of  the  processor. 

The  software  can  be  broken  down  into  two  categories: 

a.  Control  Software  -  This  software  supplies  the  linkage 
between  the  host  computer  operating  system  and  the  array  processor. 

It  is  usually  in  the  form  of  a  device  driver. 

b.  User  Software  -  This  software  enables  a  user  to  write 


programs  for  an  array  processor.  Typically,  this  can  be  done  at  two 
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levels.  A  user  can  program  in  Fortran  by  writing  a  program  for  the 
host  computer  which  consists  of  a  series  of  calls  to  the  array  processor 
math  library.  This  math  library  is  supplied  by  the  vendor  and  consists 
of  a  library  of  mathematical  routines  which  have  been  coded  for  the 
array  processor.  Figure  D  2  shows  an  example  of  such  a  program  to 
compute  a  Fast  Fourier  Transform.  Obviously,  when  the  array  processor 
is  used  in  this  manner,  its  internal  structure  is  transparent  to  the 
user. 

The  second  level  of  programming  is  to  program  the  array  processor 
directly  in  assembly  language.  Usually  the  vendor  supplies  an  assem¬ 
bler,  simulator  and  debug  aids  to  assist  the  programmer.  Figure  D.3 

shows  such  a  program  written  in  the  assembly  language  for  the  FPS 

2  2 

AR-120B.  This  program  calculates  Ci  -  Ai  +  Bi  where  i  ranges  from 
1  to  N.  The  y  axis  of  the  figure  represents  machine  cycles,  while  the 
x  axis  represents  flow  through  the  pipelines.  The  program  reduces  to 
a  4-cycle  loop.  However,  this  loop  does  demonstrate  the  parallel  struc¬ 
ture  of  array  processors,  for  example,  on  the  first  cycle  of  loop,  a 
memory  fetch,  a  floating  multiply,  a  memory  save  and  floating  add  are 
all  in  progress  on  the  same  machine  cycle  (contrast  this  with  a  conven¬ 
tional  computer).  At  this  level,  the  programmer  has  to  be  aware  of  the 
internal  structure  of  the  array  processor  to  maximize  performance. 


CALL  APCLR 


Clear  array  processor 

CALL  APPUT  Transfer  data  to  array  processor 

CALL  CFFT  Perform  complex  FFT 

CALL  APGET  Transfer  results  to  front-end  computer 

Figure  D.2.  AP  Fortran 


FETCH  STAGE 

MULTIPLY  STAGE 

ADD  STAGE 

FETCH  A 

FETCH  B 

NOP 

SAVEX 

FETCH  A; 

FMUL  A,  A,  SAVEY  B 

FETCH  B; 

FMUL  B,  B 

NOP; 

FMUL 

SAVEX  A; 

FMUL;  SAVEY  A2 

FETCH  A; 

FMUL  A,  A,  SAVEY  B 

FADD  B,  A 

FETCH  B; 

FMUL  B,  B 

FADD 

NOP; 

FMUL; 

DEC  N 

SAVEX  A; 

FMUL ;  SAVEY  A2 

STORE  C;  BGT  LOOP 
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Table  E.l.  Main  CFAR  Simulation  Program 


PBOGPAN  TC  '’"LLATE  A  CELL-AUERACING 

CPF»T£-  r  CFAP  FOP  CONPARISON  OF  THE  TD0 
ALm'PPhnS  £.  -  NTE  CARLO  TECHNIQUES •  THE 
,■  ,1NFAPIS0N  C>»'<  EE  "ODE  TOP  UARIOUS 
EN"!P0Hf1ENTs  HNO  u0*0  LENGTHS. 


SEEB  -  UN  IT OPR  RANDON  NO.  GENERATOR  SEED 
STDU  -  STANDARD  DEVIATION  Of  OAUff***  HOISE 
NEAN  -  NEAN  UAUJE  OP  GAUSSIAN  NOISE 
SNRI  -  INPUT  SIGNAL-TO-NOISE  RATIO 

NDET  -  DETECTOR0 LAUi^O.  LA“:!iKIHEW>'1 

3&  :  35  SS5l?  tfV  FALSE  ALARN 

PFDGO  -  DESIGN  GO  PROBABILITY  OF  FALSE  ALARN 

HHCR  -  WJHBER  OF  NONTE  CARLO  RUNS 

VpOF  -  RUN  INDEX)  PFA**,PD>1 

NSN  -  NO.  OF  SNR  RUNS 

NC  -  NO.  Of  CELLS  COVERED 

TAU  -  RATIO  COOEBED'NON-COUERED 

1SKP  -  CELLS  SKIPPED  BY  CFAR 

10  -  8UANTI2ATI0N)  YES>» 

HI IT  -  NO.  OF  BITS 

IUEI  -  UEIBULL  CLUTTER)  N0-8,VESM 

A  -  UEIBULL  PARANETER 

CPOU  -  UEIBULL  CLUTTER  POUER 

jDCC  -  DISK  INDEX)  CA-#,G0*1 

CALCULATED  INPUTS 

CKCA  -  CA  THRESHOLD  CONSTANT  s 

CKGO  -  GO  THRESHOLD  CONSTANT 
ANPI  -  INPUT  ANPLITUDE  AT  IF 

OUTPUTS 


pfca  -  CA  PROBABILITY  OF  FALSE  ALARN  OBTAINED 

PFCO  -  GO  PROBABILITY  OF  FALSE  ALARN  OBTAINED 

PBCA  -  CA  PROBABILITY  OF  DETECTION  OBTAINED 

roco  -  GO  PROBABILITY  OF  DETECTION  OBTAINED 

OINENSION  PFDCA(9>,PFDC0(9>*NUD2(9>,NNCR(9) 

REAL  NEAN.NNCR 

DATA  PF DCA/3*i.E -3, 3* 1 .E-4,3* I. E-5'_ 

DATA  PFPG0'8.8773E-8,8.S8®9E-2»8.318lE-2, 
B.2227E-2.8. 1B43E-2.8.5318E-3. 

• . S808E - 3 . B . 2297E - 3 ,B . 9*#8t - J ' 

DATA  NUD2''4.8.16»4,8,I«.4.8.1*y 
DATA  NNCR/3*l.BCf.3*l1J«9.3»I.EiB/ 

CALL  ASSIGNO.  'DKUCFAR.PLT\8.'N£W  > 

DEFINE  FILE  3(2*8, 128, U.JJ> 

SEE0-B.2SIB48394S 

rrou.i. 

NEAN* 8. 

SNRI*#. 

NTAR'l 


( 


Table  E.l  (cont'd) 


1 


N 


ND£T»# 

IPDF-# 

TAg.i. 

TAUSR-SORT(TAU) 

ISICP-1 

lQ-« 

NBIT-8 

IUEI-# 

A-l. 

CPOU-1. 

I DCC-8 
ICI-1 

IF (NDET .EG.# )  Q- 10X82/2. XXNB IT 
IF(NDET.EQ.l )  Q-I0/2.8XHBIT 
CALL  APCLR 

CALL  UCLR<*. 1,32767) 
r  ft[  1  ftP  UR 

CALL  APPUT(A.13,1,2> 

CALL  APPUTICPOU, 14,1.2) 

CALL  APPUT(SEED,16.1.2> 

CALL  APPUT ( STDU.S, 1,2 ) 

CALL  APPUT<NEAN,6.1,2> 

CALL  APPUT(TAUSS,ll.t,2> 

CALL  APPUT(Q, 12.1.2) 

CALL  APUD 

IFdPDF.EQ.l)  URITE(6.1#2> 

IK  F0RBAT(/1H1.2H  N,6X,  4HSNRI .  13X.SHPFDCA.  1 3X.4HPDCA, 

1  13X,4HPDGO> 

IFdPDF.EO.#)  URITE<6,1#3> 

1#3  FORflAT(  /1H1.2H  N,6X.SHPFDCA, 13X, 4HPFCA. 

1  13X.SHPFDG0,13X,4HPFG0> 

IFdPDF.EO. a>  URITE<6,1#S) 

1#5  F0RAAT(/1H1.2H  H,  3X,2HNC,6X,5HPFDCA,  13X.SHPFTCA,  13X.4MPFCA, 
1  13X.5HPFDG0, 13X.4HPFG0 ) 

C  AAIN  LOOP 

DO  E  J-3.S.3 
SNRI-#. 

NU-NUD2<J) 

Nue-Nusa 

HCLU-# 

IFdPDF.EO. 2)  NCLU-NU2M4 
CKCA-<PFDCA<J>XX<-1./NU2>-1. ) 

CXCO-(PFDCO«J)XX<-l./NU>-l.  > 

IFCNDET .EG.  t )  CKCA-S0RT<CICCA8NU2>T1 . 14/NU2 
IFINOET.EQ.l )  CKG0*S0RT(CKG08NU )t 1 .14/MU 
CALL  APFUTCCKCA, 1.1,8) 

CALL  APPUT (CX CO, 8, 1,2) 

CALL  APUR 
FNACR-NACJM  J) 

IFdPDF.EO. 1)  FNACR-1.E5 

MAON- IATIFNACR/ ( 1 .E5 AFLOAT dSKP  >- 1 .  >  > 

NSH-21 

IFCIPDF .NE. 1 )  NSN-1 
C  SIONAL-TO-NOISE  LOOP 

DO  3  I-l.NSN 

AAP I  •  8TD01  SORT  ( 8 .  X 1  • .  ( t  ( SNR  I  / 1  • .  ) ) 

CALL  APPUTl AAPI ,10,1,8) 

CALL  APWR 

NUl-NkHl 

NM14-NU1+14 

AC-0 

C  NONSTATIONARV  CLUTTER  LOOP 

DO  S  ICLU-l.NCLU 
IF < ICLU.GT . NU1  )00  TO  11 


i 
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Table  E.l  (cont'd) 


PFTCP*PFDCA'  J  )**< FLOAT (NU2-NC  >/FL0AT<NW2>  ) 

PFTCA»PFTCA*(  1  .*TPW<PFDCAI  J  >»*< -1 ./NU2  >~1 .  !  >M!  -NC  > 

CO  TO  99 

11  IF(ICLU.LT.NU14)  GO  TO  5 

PF TCP • PF DC A ( J  )  * * < F LOP T f NC  ) / F LOOT f NU2 > ) 
PFTCA*PFTCAt(l.*(PFDCA<J >**<- 1./NU2 )-l. )/TPU >**<NC-NU2 ) 

99  CONTINUE 

CPU  PPPUT(SEED.4.1,,2) 

CPU  PPPUT  ( SEED  ,7,1.2) 

CPU  APUR 
CPU  UCLR<9.1.4) 

CPU  TWIN ( NU . NNON , NTPR . NDET , ICLU , I SKP . IQ . 1 UE I ) 

CALL  apur 

CALL  APGET(DUN.9,4,5 ) 

CALL  APUD 

IFCNCUI.NE.9)  CO  TO  19 
C  PFP  AMO  PD  CALCULATION 

PFCA-CA/'(Ff94CR) 

PF CO-CO/ <FNNCR> 

PDCAII )-<DCA-CA>/FNNCR 
PDCO< I) • < DGO-GO ) /FNNCR 

IF(IPDF.EQ.l)  UR1TEI6.191 >NU2. SNRI ,PFDCA( J > ,PDCPt I I.PDGOC 1 > 
IF <  IPDF.F.Q.9 )  WRITE ( 6, 191  )NU2,PFDCA(  J  ) ,PFCP.PFDCO(  J  ) .PFGO 
191  FORNPT(2X, 12, 3XE12 .5, 4< SX.E12 .5 >) 

GO  TO  5 

C  NONSTATIONARV  CLUTTER  PFA  AND  PD  CAL. 

19  NCI’NCM 

PDCACNC1 l^CAxPHMCR 
PDCA(NC1*NU2*1  )-CO/FNBCR 
NC2-NC1*<NU2«1>*2 
PDCA(NC2)-PFTCA 

URITEC6, 194 >  NU2,NC,PFDCP(J  >,PFTCP,PDCP(NC1 >,PFDCO( J ) 

1  ,PDCP<NC1«NU2M) 

194  FORN9T<2X,ia.3X.I2.3X,E12.S.4(SX,EI2.5>> 

NC-NC+l 

5  CONTINUE 

SHRI-SNRIU. 

3  CONTINUE 

C  DISK.  WRITE 

IFUP9F.E0.9)  GO  TO  2 
IF(IPDF.EQ,2)  NSM>3S(MUa«l ) 

IF(  IDCC.E0.9)  URITEO'ICI)  IPDCPC  IP  >.  IP-l.NSN > 

IF(  IDCC.EO.  1 )  URITEO'ICI)  (PDGOC  IP),  IP-l.NSN ) 
ici-icm 
2  CONTINUE 

STOP 
END 
t 
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no 

SUBROUTINE  THLN 

PURPOSE:  To  generate  target  noise  and  clutter  inputs  to  CFAR 

program. 

FORTRAN  CALL:  Call  THLN  (NWD2,  NMON,  NTAR,  NDET,  ICLU,  ISKP,  IQ) 

PARAMETERS:  NWD2  =  Half  of  CFAR  window  width 

NMON  =  Number  of  program  passes 
NTAR  =  Target  model  number 

NDET  =  Detector  law:  Square  Law  =  0;  Linear  =  1 
ICLU  =  Number  of  cells  covered 
ISKP  =  Cells  skipped  by  CFAR 
IQ  =  Quantization:  No  =  0;  Yes  =  1 
IWEI  =  Weibull  Clutter:  No  =  0,  Yes  =  1 
EXTERNALS:  VSQRT  ,  VRAND,  VLN ,  VFILL  ,  QUANT,  RAN DM,  VSQ,  VMVL, 

VADD,  CFR,  VNEG,  WEIBULL 

SCRATCH:  SP  (0-6,  12-14),  DPX  (-4,  3),  DPY  (0,  1) 


Ill 


Table  E.Z.  Subroutine  THLN 


AP  PP06PAN  T0  GENERATE  TARGET  an: 

hC  I  St  Input  fop  ofar  program 

THERE  APE  EIGH*  INPUTS i 

MUDS  -  HALF  OF  CFAP  WINDOW  UIDTH 
MHOM  -  HO.  OF  PROGRAM  RUNS 
HTAP  -  TARGET  MODEL  MO. 

HDET  -  DETECTOR  LAUiSO.  LAW**; IJNEAR-1 
ICLU  -  MO.  OF  CELLS  COVERED 
ISFP  -  CELLS  SKIPPED  »V  CFAR 
10  -  QUANT  2  ZAT  2  OH  1  N0-*,VES*1 

IUEI  -  UEItUL:  HO*».vES*l 

•ENTRV  THLN.8 
•EXT  USIN.UCOS 

•EXT  USQRT . URAHD, VLN, UF I LL.OUAMT, UNUl 
•EXT  RAMON, USQ,USMUL,UADD.CFAR,VNEC 
MUDS  «EOU  • 

NHON  <EOU  i 
HTAR  KQU  2 
NDET  KOU  3 
ICLU  KOU  4 
ISKP  KOU  S 
IQ  KQU  6 
IUEI  KOU  7 
ICLA  KOU  13 
NlM  KOU  14 
1CLT  KOU  12 

THIMiLDDPAi  D»*13.  _ 

NOV  NTAR.HTAR;  DPXtiKSPFN 
MOV  NUD2.MUD2jDPV<*XSPFN 
NOV  NMON ,  NMOM  j  DPX  ( 1  KSPFN 
MOV  NDET.NDCT,DPX<2XSPFN 
NOV  ICLU.ICLU;DPX(-1 >CSPEH 
MOV  ISKP,  ISKPjDPXf  -2  KSPFN 
MOV  I0,IQ|DPX(3XSPFN 
NOV  IUE1,IUEI|DPYC3XSPFN 
LDSPI  14,01-10. 

NOV  MD*,Nl**iDPY<  1  XSPFN 

UTSPI  13,01*493,  _ 

SUD  NU02.ICLA,DPX(~3XSPFN 
NOV  NN0N.NN0N,DPX(-4XSPFN 

LOOP i  LOOP*!  01*13.  _ 

•  QUASSIAN  RANDOM  HO. GENERATOR 

LDSPI  S,D»-4. 

LDSPI  1 ,D1*4*D. 

LDSPI  2.0*-2*4**. 

LDSPI  3,Di*G. 

LDSPI  4.DS-S. 

LDSPI  5. Di-21***. 

LDOPAj  Dl-4. 

JSR  RAMON 
LDOPAj  01*13. 

LDSPI  7,D1*DPY(3) 

NOV  IUEI, IUEI 


Table  E.2  (cont'd) 


LDDt*.  DB-4 
JSP  .'fcasj 
LOW  n .  DB*  I  - 
LDiPi  P.Dt-2!40«. 
LDSK  1.  DB-1. 
LDSPI  2.DB-214B0. 
LDSPI  3.DB-1. 

LDSPI  4.DB-1820B. 
LDDPP;  DI-4. 

JS»  ULN 
LDDPP;  DB-13. 

LDSPI  B.DB-2148B. 
LDSPI  1 . DB-1 . 

LDSPI  2.DB-214B8. 
LDSPI  3. DB-1. 

LDSPI  *  DB-1B2B8. 
LDDPPj  L!'<. 

JSP  USQRT 
LDDPA;  DB-13. 

LDSPI  B.DB-4. 

LDSPI  I.DB-48B. 
LDSPI  2. DB-1. 

LDSPI  3.DB-1B28B. 
LDDPPj  DB-4. 

JSP  URAHD 
LDDPPj  DB • 1 3 . 

LDSPI  B.DB-4BP. 
LDSPI  l.OB-1. 

LDSPI  2.DB-1S. 
LDSPI  3.DB-4BP. 
LDSPI  4.DB-I. 

LDSPI  S.DB-1BZBB. 
LDDPPj  DB-4. 

JSP  USNUL 
LDDPPj  DB-13. 

LDSPI  B.DB-4B*. 
LDSPI  l.DB*I. 
LDSPI  B,DB-1D6BB. 
LDSPI  3, DB-1. 
LDSPI  4.DB-IKN. 
LDDPPj  DB-4. 

JSP  VSIN 
LDDPPj  DB-13. 
LDSPI  S.DB-4SS. 
LDSPI  l.DB-1. 
LDSPI  2.DB-4SS. 
LDSPI  3 . DB • 1 . 
LDSPI  4.DB-1B2BP. 
LDDPPj  DB-4. 

JSP  UCOS 
LDDPPj  DB-13. 
LDSPI  S.DB-4BS. 
LDSPI  l.DB-l. 
LDSPI  2.DB-214BB. 
LDSPI  3.DB-1 . 
LDSPI  4.DB-4SB. 
LDSPI  S.DB-I. 
LDSPI  6,DB-1B2BB. 
LDDPPj  DB-4. 

JSP  VMUl 
LDDPPj  DB-13. 
LDSPI  B.DB-ISSBS. 
LDSPI  1,DB-1. 
inai  ».Dft«ai4tiP. 


Table  E.2  (cont'd) 


LDSPI  3. Cl*  1  . 

LDSPI  4.DI-10600. 

LDSP !  5,01*1 . 

lds°i  6,Di-ie2«*. 

LDDPA;  01*4. 

JSR  UHUL 

LNWE 1 1  LDDPA;  Di* 1 3 

LDSPI  8.DB-DPXI0J 
MM  NTAR.NTAR 
BCD  CASS 
JNP  CAS1 

•  CASE  •  TARGET  GENERATOR 
CAS#!  LDDPA)  DB-13. 

LDSPI  0.DI- 10  - 
LDSPI  1, DB-30000. 

LDSPI  8.DB-1. 

LDSPI  3,DB*!M#e. 

LDDPAi  Dl*4. 

JSR  WEILL 
LDDPAi  DB*13. 

JAP  QCH 

•  CASE  1  TARGET  GENERATOR 
CASH  LDDPAi  DB-13. 

LDSPI  0.DB-7. 

LDSPI  1 , DB ■ 30000 . 

LDSPI  2.DB-1. 

LDSPI  3,DB-1000B- 
LDDPAi  DB*4 . 

JSR  URAND 
LDDPA;  DB-13. 

LDSPI  #.DB-30##B. 

LDSPI  1 .DB* 1 . 

LDSPI  2.DB-30000. 

LDSPI  3.DB-1. 

LDSPI  4.DB-10000. 

LDDPAi  DB-4. 

JSR  WIN 
LDDPA)  DB-13. 

LDSPI  0, DB-30000. 

LDSPI  I.DB-1. 

LDSPI  S.DB-30000. 

LDSPI  3.DB-1. 

LDSPI  4.DB-10000. 

LDDPAi  DB-4. 

JSR  WNEC 
LDDPA)  DB-13. 

LD6PI  0.DB-30000. 

LDSPI  l.DI-l. 

LDSPI  8. DB-30000. 

LDSPI  3.DB-1. 

LDSPI  4, DB- 10000. 

LDDPAi  DB-4. 

JSR  WSORT 
LDDPA)  DB-13. 

LDSPI  0, DB-30000. 

LDSPI  l.DB-l. 

LDSPI  8.DB-10. 

LDSPI  3, DB-30000. 

LDSPI  4.DI-1. 

LDSPI  D.DB-10000. 

LDDPAi  DB-4. 

JSR  WDNUL 
0CHI  LDDPAi  DB-13. 

•  NOISE  ♦  TARGET 
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Table  E.2  (cont'd) 


ldspi  a.  oioaata. 

LDSPI  l.OB-l. 

ldspi  e,c»-sa«. 

LDSPI  3.01*1. 

LDSPI  4.0B-3MBB. 

LDSPt  5,00*1 . 

LDSPI  E.OBMMae. 

LDOPAi  01*4. 

JSB  UPDO 
LDDPAj  DO* 13. 
idspi  a.DB*3Maa. 

LDSP1  l.DB-1. 

LDSPt  e,o**3«aae. 

IDSPI  3.01*1. 

LDSPt  4,DB-1««M. 

LDDPftj  01*4. 

JSB  WSQ 

•  NOISE  ONLV 

LDOPPj  DB* 1 3 . 

LDSPt  «,DB-DPX(B) 

LDSPI  13.DB-493. 

SUB  NUD2,  ICLPiDPX  t  -3  XSPFN 
LDSPI  12,DB-DPX<-1) 

CLUt  LDDPBjDB*13. 

LDSPI  B,D8*DPX(-3> 

LDSPt  l,DB*DPX(-2) 

LDSPI  2 , DB • 1 1 . 

LDSPI  3, DB-DPXC -3 ) 

LDSPI  4, DB-DPXf -2 ) 

LDSPI  S.DB-2M. 

LDDPPi  DB-4. 

JSB  USHUL 
LODPBj  DI-13. 

LDSPI  13, DB*DPX(-3 ) 

INC  ICLP|DPX(  -3  XSPFN 
DEC  ICLT 
BED  NOIS 
JNP  CLU 

NOISx  LDDPB,  DB*13. 

ldspi  a,DB*4ea. 

LDSPI  t , DB • 1 . 

LDSPI  2.DB-4M. 

LDSPI  3,DB*1. 

LDSPI  4.DB-2B4SS. 

LDDPPi  DB*4. 

JSB  VSO 
LOOP*,  DB-13. 

LDSPI  S,DB-3MM. 

LDSPI  I.DB-1. 

LDSPI  S.DB-1B7SS. 

LDSPI  3.DB-1. 

LDSPI  4.DB-3MSS. 

LDSPI  S.DB-t. 

LDSPI  S.DB-1MM. 

LDDPftj  DB-4. 

JSB  UBDD 
LDDPBj  DB-13. 

LDSPI  S.DB-4M. 

LDSPI  l.DB-1. 

LDSPI  2,DB-1KM. 

LDSPI  3.DB-1. 

LDSPI  4.DB-4SS. 

LDSPI  S.DB-l. 

LDSPI  S.DB-1B2DB. 


Table  F.2  (cont'd) 


LDDPA,  DB*«. 

J SR  vftDD 
DETECTOR  LAw 
LDDPA,  DB • 1 3 . 
LDSPI  3.DB-DPXI8; 

nou  ndet.mdet 
BMC  END 


JNP  SLD 
•  LINEAR 

END'  LDDPA:  DB 
LDSPI  0 


13. 

DB*3M00. 
LDSPI  l.DB-1. 

LDSPI  2,DB*3MBB. 
LDSPI  3.DB-1. 

LDSPI  4.DI-1MM. 
LDDPA,  DB*4 . 

JSR  USQRT 
LDDPA;  DB* 13 . 

LDSPI  »,DB*4e». 
LDSPI  l.DB-1. 

LDSPI  S,DB*4«e. 
LDSPI  3, DB* 1 . 

LDSPI  4.DB-102S*. 
LDDPA.  DB*4. 

JSR  USQRT 
SQUARE  LAV 
SU»  LDDPA  j  DB*13 . 

QUANTIZATION 
LDSPI  6.DB-DPX13) 
NOU  IQ, 10 
BNE  LOOPS 
JAP  L00P1 

LOOPS:  LDDPA,  DB-13. 

LDSPI  0.DB-40O. 
LDSPI  1.DB-1S. 

LDSPI  S,DB*10SM. 
LDDPA, DB*4. 

JSR  QUANT 
LDDPA,  DB*13. 

LDSPI  0,DB*3Q000. 
LDSPI  1,DB*IS. 

LDSPI  2,DB*10Q00. 
LDDPA, DB* 4. 

JSR  QUANT 
NOISE  ONLV  CFAR 
LOOP1 I  LDDPA,  Dl • i 3 . 

LDSPI  0,DB*DPYt0> 
LDSPI  1,DD*1SM0. 
LDSPI  2.DB-S. 

LDSPI  6,08*493. 
LDSPI  7,DB*S97. 
LDSPI  1S,DB*5S0. 
LDSPI  13,DB*DPX(-2) 
LDDPA,  DB*4. 

JSR  CFAR 

TARGET  *  NOISE  CFAR 
LDDPA,  DB*13. 

LDSPI  t,DB*DPVIS  > 
LDSPI  1,DB*1MM. 
LDSPI  2,DB*S. 

LDSPI  6,DB*493. 
LDSPI  7,DI*6S7. 
LDSPI  1S.DB-3S000. 
LDSPI  I3,DI*DPX(-S> 


SUBROUTINE  CFAR 


PURPOSE: 

FORTRAN  CALLS: 
PARAMETERS: 

EXTERNALS: 


I 


To  simulate  two  CFAR  algorithms:  Cell  averaging 
and  greatest-of 

Call  CFAR  (NWD1,  NT,  OA,  NN,  TB,  TA,  TCI) 

NWD1  =  Half  of  CFAR  window  width 

NT  =  Total  number  of  inputs 

OA  =  AP  output  address 

NN  -  Number  of  cells  skipped 

TB  =  AP  address  before  cell  of  interest 

TA  =  AP  address  after  cell  of  interest 

TCI  =  AP  address  for  cell  of  interest 

None 

SP  (0-2,  6-7,  13-15) 


SCRATCH: 


Table  E.3.  Subroutine  CFAR 


•TITLE  CFAR 

•  AP  PROCRAP  TO  SIMULATE  TUO  CFAR 

•  ALGORITHMS)  CELL  AUCRACING  AND 

•  GRCATEST-OF.  THERE  ARE  SEUEH 

•  INPUTS: 

•  MUDS  -  HALF  OF  CFAR  WINDOW  WIDTH 

•  NT  -  TOTAL  NO.  OF  INPUTS 

•  OA  -  AP  OUTPUT  ADDRESS 

•  nn  -  NO.  OF  CELLS  SKIPPED 

AP  ADDRESS  REFORE  CELL  OF  INTEREST 
AP  ADDRESS  AFTER  CELL  OF  INTEREST 
AP  ADDRESS  FOR  CELL  OF  INTEREST 


•  T» 

■  TA 

•  TCI 
•ENTRV  CFAR,? 

NUDS  SEOU 
NT  SEOU 


1 

oa  seou  a 

NN  SEOU  13 

TB  SEOU  6 

TA  SEOU  ? 

TCI  SEOU  IS 

NS  SEOU  M 

CFAR!  NOW  NUDS, NS 
DEC  NS 
SUB  NS, TB 
LDTNAj  DB* 'ONE 
NOP 

DPXC2XTN 

NOU  OA, OA j  SETNA 


U 

'A 

< 

i 

J 

■  •< 


DPVO  )<ND|  INCNA 

NOP 

HOP 

DPV(-4XHD 
LDSPI  1S.DB-B. 
NOU  IS, IS)  SETNA 


DPY(SXND)  INCNA 


DPV(-3XHD 

•  LOOP  FOR  AMINO  N  PREVIOUS  CELLS 
LOOPS*  FAM  ZERO, ZERO | NOU  TB.TBj SETNA 

FAM 

NOU  NUD8.NS 
LOOP1 *  INCNA;  DPXCND 

FAM  DPX.FAiDEC  NS 
FAM  )BOT  LOOP1 
DPX(-EXFA 

•  LOOP  FOR  ADDING  N  FORWARD  CELLS 

FAM  ZERO.ZEROiNOU  TA.TAj SETNA 
FAM 

NOW  NUD2.NS 
LOOPS*  INCNAiDPXCHO 

FAM  DPX.FAtDEC  NS 
FAM  i  DOT  LOOPS 
DPVUXFA 

•  ADDING  SUNS  AM  STORING  IN  DPX(S) 

FAOD  DPVUI.DPXC-S) 

FADD 

DPXKFA 

•  FINDING  LARGEST  SUN  AMD  STORING  IN  DPX<-3> 


Table  E.3  (oont'd) 


fsvjs  DPX(-e>,DPvcn 
FSUS 

0Px<-3XDP*'-a) 

BFGT  LOOPS 
DPX(  -3  XDPYS  1 ) 

LOOPS l  ADD  MM, TA 
ADD  Mh.Tl 

•  MULTIPLY  SUMS  *Y  X/N 
FMUL  DPX(S),DPY<0) 
FNUL  DPX< -3 ), DPY t -3 ) 
rnuL 

FNULjDPXCFN 
FNULt  DPX(-3XFn 
nou  TCI.TCIjSETMA 
NOP 


NOP 

DPY(2XND 

«  PERFORM  CA  THRESHOLD  COMPARE 
FSUI  DPY( 2 )»DPX(#> 

FADD 

NOP 

BFCE  CROSS 
SR  FIN 

CROSS >  FADS  DPX(2>.DPY(3> 

FADD 

DPVOKFA 

*  PERFORM  CO  THRESHOLD  COMPARE 
FINS  FSUD  DPY(2>,DPX(-3> 

FADD 

NOP 

SFGE  CROS1 
SR  FIN1 

CROSlt  FADD  DPXC2>.DPV(-4> 

FADD 

OPYC-4  XFA 
FINIS  ADD  NN.TCI 
SUB  NN.NT 

sea  fin2 

JNP  LOOP'S 

•  STORE  NO.  OF  DETECTIONS 

FIN2S  NOU  OA.OA 

NICDPV(3)|SETNA 

NI<DPY< -4 )j INCNA 
NOP 

LDDPAi  DS-4. 

RETURN 


SEND 

S 
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PURPOSE: 

FORTRAN  CALL: 
PARAMETERS: 


FORMULA: 


EXTERNALS: 

SCRATCH: 

NOTES: 


SUBROUTINE  VRANDX 

To  generate  an  array  of  random  numbers  uniformly 
distributed  between  0  and  1. 

Call  VRANDX  (A,X,I,N) 

A  =  Address  of  starting  seed 

X  =  Base  address  of  output  array 

I  =  Increment  between  elements  of  output  array 

N  =  Number  of  output  samples  desired 

Technique  used  is  multiplicative  congruential  method. 

X(0)  =  M0D(8*A,1 -0)  where  B  =  27.0 

X(M)  =  M0D(B*X(M-1),  1.0)  for  M=1,2...,N-1 

None 

SP  (0-3,  DPX  (0-2),  DPY  (0) 

1.  Preferred  starting  seed  is  0.2510637948. 

2.  At  completion  the  seed  is  set  for  the  last 
number  generated. 
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Table  E.4.  Subroutine  VRANDX 


URANDX  -  VEGTOP  RANDOM  NUHBERS  /CC"NON'  •  RE’-  3.9.  AUG  77  SSttS* 

•  FOR  EITHER  MEMORY 
•TITLE  UP AND* 

•ENTRY  URANDX. 4 

• 

•  - ASSTffACT _ 

•FILLS  UCCTOR  C  WITH  A  SEOUENCE  OF  FLOATING  POINT  RANDOM  NUMBERS 
•UNIFORMLY  DISTRIBUTED  BETWEEN  AND  l.B.  SEOUENCE  IS  GENERATED 
•USING  A  SEED  A.  FOLLOWING  GENERATION  THE  SEED  IS  SET  TO  THE  LAST 
•RANDOM  NUMBER  GENERATED.  THUS  ALLOWING  THE  SEOUENCE  TO  IE  CONTINUED 
•IN  THE  NEXT  CALL  TO  URAMP.  SUGGESTED  SEED  FOR  FIRST  CALL  IS  B.2S1B637948. 

t 

•TECHNIQUE  USED  IS  MULTIPUCATIUE  CONGRUENT  I AL  METHOD. 

’FORMULAS  C C • > «MOD ( DBA .  1 .  • )  WHERE  B-27.0 

•  C(MK  )*HOD( BSC< (H-l)Kl.l.Q)  FOR  M*1  TO  N-l 

•  AND  A  •  C((N-nm 

• 

•  STATISTICS-  — 

•LANGUAGES  ARAL 

•EQUIPMENT)  AP-12B  WITH  CITHER  MEMORY 
•SIZES  IS  LOCATIONS 
•SPEEDS  INTROS  S  CYCLES 

•  LOOPS  7-8  CYCLES  <7.1  CYCLES  AVERAGE) 

•  COLUMNS /LOOP)  1 

•  FLOPS/LOOPl  3 

•  1.19N  ♦  8.83  USEC,  FOR  167  NSEC  CLOCK 

•  NECAFLOPSs  2.52 
•SUMOUTINES  USED:  NONE 
•AUTHORS  R.S.  NORIN 
•SATES  JAN  77 

• 

•  USAGE 

•FORTRAN I  CALL  URANDXt A.C.K.N) 

•APAL>  JSR  URANDX 

■S-PAO  PARANETERS 

•  NAME  NUNBER 

A  8COU  •  'ADDRESS  OF  SEED 

C  SEQU  1  'BASE  ADDRESS  OF  DESTINATION  UECTOR  C 

K  SEQU  2  • INCREMENT  BETWEEN  ELEMENTS  OF  C 

N  SEQU  3  ‘NUMBER  OF  ELEMENTS  IN  C 

•TABLE  MEMMV 
ONE  8EBU  'ONE 

• 

•SCRATCMs  SP(I,3>,  DPX(B-2).  DPYCB),  DPA  UNCHANGED 


URANDX I  MOW  A.Ai  SETNA 
RPSF  B>  DPX<DB 
LDTNAiDB-ONE 
RPSF  FMASKt  DPX(2KDB 
FNUL  DPX.itti 
DPXtl  KTMj 
SUB  K.C 

LOOP*  FNUL 
FNUL 

FSUB  FN.DPXU  J, 


DPVCFM 

FADD 

FAMD  DPX  <  2 ) , DPV 
FADD  ZERO. DPV i 


•GET  SEED  A 
•GET  MULTIPLIER  B 

•GET  FRACTION  MASK 
■•BA 

•SAVE  !.• 

•BACK  UP  DESTINATION  ADDRESS 

•PUSH 

•PUSH 

•FORM  ISA-1  SINCE  NETHOD  OF 

•EXTRACTING  FRACTION  WILL  BE  DIFFERENT 

•IF  ISA<1.B 

•SAVE  DBA  FOR  LATER 

•PUSH 

•ASSURE  ISAM  SO 

•FRACTION  CAN  BE  EXTRACTED  WITH  MASK 
•OCT  FRACTION  DIRECTLY  IF  ISA<1 


i 
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Table  E.4  (cont'd) 


DtC  N j 
BFGT  K 


iTl 


FA DDi  MOU  N.N 

ctii  add  k.c.setmaiHKfaj 

FNUl  DPX.FAj 
BNE  LOOP 

DOME!  NOLI  A. A)  SETMA  j  MI <FAj 


RETURN 

B  i  VP  27. • 

FflASK:  VP  • .  999999993S 
•EMC 

* 


•DECREMENT  COUNT 

*CET  FRACTION  IMMEDIATELY  FROM 

•Fa  IF  |tA>l 

•This  IS  an  EXTRA  CYCLE 

•In  LOOP  If  |*a<i 

•STORE  RANDOM  NUMBER 

•START  FORMING  NEXT  NUMBER 

■CONTINUE  UNTIL 

•DONE.  THEN  STORE  LAST  RANDOM 

•NUMBER  AS  The  NEW  SEED. 

■Then  EXIT. 

•MULTIPLIER  CONSTANT 
•fraction  MASK 
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PURPOSE: 

FORTRAN  CALL: 
PARAMETERS: 


FORMULA: 


EXTERNALS: 
SCRATCH: 
NOTES : 


SUBROUTINE  RAN DM 

Generates  an  array  of  random  numbers  which  are 
independent  and  have  a  Gaussian  distribution. 

Call  RAN DM  (A,  X,  N,  MEAN,  STD,  SCR) 

A  -  Address  of  starting  seed 
X  =  Base  address  of  output  array 
N  =  Number  of  samples  desired 
MEAN  =  Location  of  the  desired  mean  value 
STD  =  Location  of  the  desired  standard  deviation 
SCR  =  Base  address  of  scratch  storage  (N  words 
of  scratch  storage  are  required) 

Starting  from  two  random  numbers  u^  and  u ^  which  are 
uniformly  distributed  between  0  and  1,  two  Gaussian 
numbers  with  desired  mean  and  standard  deviation  are 
obtained  as 

m-|  =  n-|  a  +  p 

m2  =  n2  a  +  ^ 

where 

n-j  =  /-2  In  u-j  cos  Um^) 
n^  =  /-2  In  u.|  sin  (En^) 

VRANDX,  LN,  SQRT,  COS,  SIN 

SP  (0-9),  DPX  (-3,  -2,  0-2),  DPY  (-3,  -2,  0) 

N  words  of  scratch  storage  are  required 


Table  E.5.  Subroutine  RANDM 


•title  randb 

•£NT»y  RanDB.6 
•EXT  URAHDX.LH.SORT.CCS.S'm 
SEED  (EOJ  » 

X  SEOU  1 
N  SEQU  e 
BE  AH  SEQU  3 
STDU  SEQU  4 
Z  SEOU  8 
HS  SEOU  6 
XS  SEOU  7 
N2  SEQU  1® 

V  SEOU  11 
XI  SEQU  3 
X2  SEOU  4 

RANDNt  BOU  BE  AH, BE  AH;  SETBA 
BOU  N,HS 

BOU  STDU, STDU,  SETBA 

INC  H;  DPX !  -3  XBD 

NOUR  N,N2,  LCTBA;  D*-'TUOPI 

BOU  X, V;  DPV< -3 ) CBD 

BOU  X.XS;  DPX(-2KTB 

ADD  N2.V 

BOU  N,3 

LDSPI  2.DS-1. 

JSR  URAHDX 
BOU  XS.XiSETBA 
BOU  XS.Xl 
DEC  XI 
BOU  N2.N 

LOOPli  DPXCBD;  JSR  LH 
JSR  SORT 
INC  X, SETBA 
DEC  N 

INC  XI,  SETBA,  BKDPX,  SHE  LOOPl 

BOU  N2.N 

BOU  V.X, SETBA 

BOU  Z.Xl 

DEC  XI 

LOOPS!  FBUL  DPX(-2),BD 
FRUL,  BOU  XI, X2 
FBUL,  ADD  NS.X2 
DPY1-2XFB 
DPXCFB,  JSR  COS 
INC  XI;  SETBA,  BKDPX 
DPX<0PY(-Z>,  JSR  SIN 
INC  X, SETBA 
DEC  N 

INC  X2;  BKDPX,  SETBA,  SHE  LOOP2 
BOU  N2.N 
DEC  V 

BOU  XS.Xi  SETBA 
BOU  XS.XS 
BOU  Z.Xlj  SETBA 
DPXCBD,  DEC  X2 
ADO  BS.X1,  SETBA 
LOOP3I  FBUL  DPX.BD 

FBUL,  INC  X,  SETBA 
FBUL  DPX.BD 

INC  X2,  SETBA,  BKFH;  FBUL 

FBUL;  DPX CBD 

INC  Vj  SETBA,  BKFH 

NOP 

INC  Z,  SETBA 


KF 


Table  E.5  (cont'd) 


DEC  Nj  LD’NA,  D8-  SOB"? 

INC  XI,  SE’Pfi;  »NE  L',OP3 

FNUL  T*l,  DPv  •  -3  > .  no<-  /S.*;  SET1W. 

fnul:  dec  xs 

FAUL;  INC  X,  SETUP 
D»»V(-3)<rn!  FNJL  FN.BD 

fnul 

LOOPS:  FAUL  DPVI-3J.AD 

FADD  FA,DPX(-3);  INC  X,  SCTAA 

FADD;  FAUL:  DEC  NS 

INC  XS;  SETAA;  AICFA,  8NE  LOOPS 

RETURN 

SEND 

f 
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PURPOSE: 


FORTRAN  CALL: 


FORMULA: 


EXTERNALS: 

SCRATCH: 

NOTES: 


SUBROUTINE  WEIBUL 

To  generate  an  array  of  clutter  amplitudes  (one  for 
each  clutter  cell  of  interest)  where  the  spatial 
statistics  of  clutter  power  are  described  by  a  Weibull 
distribution  and  zero  correlation. 

Call  WEIBUL  (S,X,B,NC,A) 

S  =  Address  of  seed  for  random  number  generation 
X  =  Base  address  of  output  array 
B  =  Base  address  of  array  [B(i)]  containing  the 
median  powers  from  each  clutter  cell 
NC  =  Number  of  samples  (clutter  cells)  desired 
A  =  Address  of  Weibull  parameter,  a 

X(i)  =yBii7[(ln  ^)]a  i=0.1 - - (NC-1 ) 

where  [u ( i ) ]  is  a  set  of  random  numbers  which  are 
uniformly  distributed  between  0  and  1. 

VRANDX ,  LN,  EXP,  VSQRT 

SP  (0-8,  13-15),  DPX  (-4-3),  DPY  (-4-3) 

The  array  [B(i)],  which  is  prestored,  must  reflect 
the  variation  of  clutter  power  between  different 
clutter  cells. 
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Table  E.6.  Subroutine  WEIBUL 


rriTLt  yei»uu 

•ENTPV  «CI»U1.5 
*CXT  OMNOX.UN.EXP 
SEED  SEOU  • 

CO H  *EOU  8 
K  K«J  1 
N  SEQU  3 
A  «£OU  4 
CONS  SEOU  5 
rs  SEQU  6 
ns  teou  7 
KSS  KOU  4 

“tI*UL’!SS  &:c«“ 

w  5:51s,  Dev(-3X«t> 
tOSPI 

JSS  VHMiVX 
-  DEC  CONS 

NOU  SS.ICSS 
NOU  KS.K)  SETNA 
DEC  KS 
BOO  NS.N 

LOOPI DPX( • KTO;  JS*  UN 

FNOLLDPX(*l.DeV(-3>l  N 00  CONS. CONS »  SETNA 

FHUL 

FNUL 

DPV( -2  XTO 

FHUL|  INC  I tj  SETNA 
l55cL4sf*TNA,  NUFN,  »NE  UOOP 
HOP 


SUBROUTINE  QUANT 


PURPOSE:  To  truncate  detector  outputs 

FORTRAN  CALL:  Call  QUANT  (X,Q,N) 

PARAMETERS:  X  =  Base  address  of  input  and  output 

Q  =  Base  address  of  LJB  level 

N  =  Number  of  inputs 

EXTERNALS:  Div 

SCRATCH:  SP  (0-4),  OPX  (-1,  0,  1,  3),  DPT  (-1,  -2) 
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Table  E.7.  Subroutine  QUANT 


•TITLE  OUANT 
•ENTRY  GOANT, 3 
•EXT  DIV 

X  SEOU  • 

0  SEOU  1 
M  SEOU  2 
V  «£QU  3 
T  SEOU  4 

OUANTjMXJ  X.V,  SETNA 
LDTNA,  Cl"  1  ONE 

MX I  Q.O:  sctna 
DEC  Vi  DPXC-1XRD 
TABS  H>xi-l)|  DPY(*XTN 
FADD.  DPVdXNDj  DPXCeXND 
DPV(-i  XFA 

jsp  diu 

FMJL  DPXCB).DPVC-1> 

LDTNA j  DB-IHALF 
FMJL j  DPV(#KDPX(-1> 

FMJL;  DPXdXTN 
LBSPI  Ti  DB-27.J  DPXC2XFN 
LOOPlFIXT  DPXC2) 

FADD 

DPXOXFAj  INC  X,  SETUP 
F ADD  ZERO. NDPXC  3  )  j  NOV  T.T 
F  ADD 

FADD  DPXdJ.FA;  DPVCDXND 
FA»S  DP VC • ) 

FMJL  DPVdJ.FAj  FSUB  DPX < - 1 ) , ZERO 
FMJL  DPXCBJ.F'Ai  FADD 
fmjLi  Dpxc-ixbpvce) 

DPV(2)<FNj  »FCT  POS 
FSUBR  APVC2),  ZERO 
FADD 

0PV(2XfA 
POS  1  FMJL  (  DEC  N 

INC  VjSETKAi  NI<DPV«2) j  INC  LOOP]  DPXC2XFN 
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